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CHAPTER  1 


INTRODUCTION 

Laser  induced  damage  in  optical  components  is  the  limiting 
factor  in  the  design  of  high  power  lasers.  This  statement  applies 
particularly  to  coated  optical  components  due  to  the  coatings'  lower 
damage  thresholds.  The  probles  has  existed  since  the  beginning  of  the 
laser  and  is  exhibited  through  tha  research  published  in  this  field.  A 
prime  example  of  this  is  the  National  Bureau  of  Standards'  Special 
Publication  of  the  Proceedings  of  the  Symposium  on  Optical  Materials 
for  High  Power  Lasers,  or  Boulder  Damage  Symposium  which  has  .been  held 
annually  for  the  past  seventeen  years  (since  1969} . 

Understanding  of  this  phenomena  has  been  limited  because  of  the 
complexity  of  the  many  processes  involved  and  the  lack  of  the  ability 
to  model  intense  field  interactions.  In  addition,  material  technology 
is  just  beginning  to  reada  a  stage  where  adequately  pure  materials  can 
be  produced  for  high  intensity  optical  applications.  The  deposition  of 
optical  coatings  lags  in  this  area  due  to  the  extreme  complexity  of 
producing  a  quality  optical  thin  film.  The  characteristics  (packing 
density,  porosity,  purity,  stoichiometry,  uniformity,  etc.)  of  optical 
coatings  seldom  approach  their  bulk  counterparts.  The  particular 
interest  in  optically  thin  films  is  because  of  their  relative  sensi¬ 
tivity  to  intense  fields  and  the  feet  that  most  optical  components  are 
coated  to  suppress  or  enhance  reflection. 


The  conditions  of  deposition  (partied  vacuum,  temperature, 
molecular  energies,  etc.)  along  with  substrate  processing  (polishing, 
ate.)  used  in  the  production  of  optical  quality  coatings  are  also 
conducive  to  film  imperfections.  Furthermore,  molecular  growth  of 
films  not  only  tends  to  allow  trapping  of  undesirable  substances,  but 
may  even  promote  it  [Bama,  Radoczi  and  Reicha  1985].  Imperfectly 
cleaned  substrates  may  also  cause  misalignments  of  growth  (nodules) 
[Guenther  1981]  with  poorly  bonded  boundaries  and  defects  of  various 
sorts.  There  is  significant  evidence  [Hurt  and  Decker,  1984]  that 
these  inperfect  regions  are  sources  of  the  initiation  of  laser  induced 
dacage  in  optically  thin  films. 

in  practice,  optical  components  are  often  coated  witivmany 
layers  of  films.  This  multilayer  system  introduces  many  new  problems 
which  are  not  related  to  the  basics  of  the  interaction  of  the  field  and 
the  film.  For  example,  the  field  distribution  within  a  stack  of  films 
must  be  computed;  the  difficulties  of  controlling  the  quality  and 
parameters  increases  significantly  with  an  increase  in  the  number  of 
layers;  furthermore,  the  many  interfaces  (a  weak  point  of  coatings)  of 
different  materials  adds  to  the  complication.  For  this  reason,  these 
studies  are  restricted  to  a  single  layer  film. 

There  are  several  good  data  bases  for  laser  induced  damage  in 
single  layer  optical  coatings  [Walker,  Guenther  and  Nielson  1981; 

Bettis  1975;  Rainer,  et  al.  1982].  Due  to  the  difficulty  of  control¬ 
ling  and  characterizing  all  of  the  parameters  (film  and  laser) 
involved,  a  large  data  base  is  desirable  to  minimize  the  effect  of 
statistical  fluctuation  in  the  damage  threshold. 


Hie  goal  of  this  research  is  the  development  of  a  model  of 
laser  induced  damage  in  optical  coatings,  based  upon  thermodynamic 
considerations,  and  the  correlation  with  of  the  observed  trends  in  the 
data.  From  the  model  a  microscopic  theory  is  then  developed  which 
produces  the  required  conditions  for  the  thermodynamic  model. 

Most  analytic  research  in  this  field  has  been  restricted  to 
laser  interaction  with  ideal  intrinsic  materials.  These  unfortunately 
bear  no  resemblance  to  real  world  optical  coatings.  The  theoretical 
projections,  when  compared  to  data,  make  this  glaringly  apparent.  Most 
other  works  concerned  with  extrinsic  processes  have  either  concentrated 
upon  mechanists  of  absorption  or  thermodynamic  treatments.  Several 
models  have  been  developed  based  upon  absorption  and  thermal  diffusion 
that  takes  place  in  the  optical  materials  IBloeobergen  1972;  Walker,  et 
ad.  19&1;  Rcmolov  1982] .  Hone,  to  this  point,  have  demonstrated  a 
relationship  between  the  films  thermal  properties  and  the  damage  that 
takes  place.  In  addition,  none  of  the  thermal  models  have  described 
mechanises  for  the  absorption  that  takes  place  other  than  that  there 
must  be  an  absorbing  impurity  of  sane  sort. 

That  is  not  to  say  that  the  ultimate  solution  of  laser  induced 
damage  in  optical  coatings  has  been  found.  A  physically  feasible 
process  has  been  described  and  justified,  however,  which  is  completely 
consistent  with  the  most  thorough  data  sets  of  laser  induced  damage  of 
optical  coatings. 

The  only  previous  relationship  found  for  laser  induced  damage 
to  optical  coatings  with  material  properties  was  that  of  Bettis :  [1975] 
scaling  with  index  of  refraction.  Though  Bettis'  model  was  simple,  it 
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|  scaled  quite  well  with  seme  balk  damage  data,  in  addition  to  sane 

coatings.  It  did  not  ,  however,  scale  well  with  tench  of  the  data  used 
in  this  present  study.  Furthermore,  it  was  based  upon  the  classical 
analog  of  the  transition  probability  intrinsic  to  the  respective  coat¬ 
ings  and  should  not  apply  well  to  the  extrinsic  damage  observed  in  the 
coatings  of  interest.  Furthermore  it  does  not  explain  the  observed 
morphology  in  the  optical  coatings  of  interest. 

It  is  clear  to  this  author  that  there  are  many  processes  of 
laser  induced  damage  in  optical  materials  that  apply  to  cany  different 
situations.  The  processes  described  within  may  be  more  or  less 
applicable  than  described.  The  calculations  demonstrate  that  they  must 
at  least  be  a  contributing  factor. 

The  morphology  of  damage  in  fluoride  and  seme  oxide  films 
dictate  that  damage.,  and  thus  anomalous  absorption,  occurs  at  isolated 
microscopic  sites  [Walker,  et  al.  1S81] .  Therefore,  a  natural  model  is 
that  of  a  localized  site  absorbing  at  many  times  the  rata  of  the 
surrounding  film.  Due  to  the  localization  of  the  damage  and  the  scales 
involved  (diameter  of  damage  site  2a  -v,  0.1  microns,  pulse  length  5  to 
40  ns)  significant  thermal  diffusion  of  energy  must  take  place  (i.e. 

<5  *  *Dtp  >  a) .  Here  D  is  the  thermal  dif  fusivity  in  cm2/s  ana  tQ  is 
the  laser  pulse  length.  Thus,  it  is  likely  that  thermal  diffusion  of 
the  absorbed  energy  is  an  important  factor  in  the  determination  of  the 
required  incident  laser  energy  which  causes  damage. 

The  morphology  of  damage  suggests  that  an  appropriate  model  to 
try  is  that  of  a  spherical  region  which  absorbs  rac  iation.  The  scales 
involved  necessitate  chat  one-  dominant,  process  in  the  localized  energy 
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deposition  leading  to  damage  is  diffusion.  This  mathematical  model  is 
developed  in  the  following  chancer.  The  requi resent  of  a  spherical 
region  is  unnecessary.  The  spherical  g  ease  try  is  primarily  used 
because  of  the  nathematicai  tractability. 

Thermal  heat  balance  demonstrates  that  for  a  given  distribution 
of  absorbing  region  sizes  there  is  one  size  of  region  most  likely  to 
damage.  The  characteristic  size  of  the  most  easily  damaged  region  is 
determined  by  the  material  parameters  and  the  pulse  length.  This 
concept  is  applied  in  a  oroad  sense  because  even  if  the  regions  evolve 
in  time  they  have  some  average  size  over  the  duration  of  the  pulse 
which  determines  the  quantity  of  energy  absorbed.  This  concept  is 
consistent  with  the  observed  morphology  [Walker,  et  al.  1981] v'  In 
particular,  the  damage  regions  observed  are  larger  at  longer  pulse 
lengths.  The  reasonable  criterion  chosen  for  damage  is  that  the 
boundary  of  the  region  reaches  some  characteristic  phase  change 
temperature.  The  melting  temperature  is  chosen  for  these  studies  due 
to  the  observation  of  a  region  bordering  the  pits  where  material  flow 
seams  to  have  occur ed. 

The  appropriate  heat  conduction  equations  are  then  formulated 
and  solved  via  Laf  Lace  Transform  [Goldenberg  and  Tranter  1952] ,  or  via 
the  more  general  integral  transfers  [Gzisik  1980 i .  The  quantity  of 
interest  is  the  incident  laser  energy  density  which  causes  the  local¬ 
ized  region  to  heat  up  to  the  damage  criterion  point.  This  quantity  is 
exhibited  in  the  icn-hcmogenecus  source  term  of  the  thermal  diffusion 
equation.  It  is  multiplied  by  an  absorption  cross  section  for  the 


absorbing  region  and  indicates  what  portion  of  the  incident  bean  is 
absorbed  and  deposited  in  the  localized  voiuae. 

The  classical  absorption  cross  section  is  computed  frctn 
electromagnetic  theory  {Mie  Theory)  [Sorn  and  Wolf  1980}  assuming  an 
highly  absorbing  localized  region  (a  sphere) .  The  size  of  the  sphere 
used  is  that  of  the  size  of  the  region  -rest  likely  to  damage  due  to 
thermal  considerations. 

With  this  model  all  parameters  of  import  may  be  specified  with 
the  exception  of  the  properties  of  the  absorbing  region.  It  is  shown 
in  the  following  chapter  that  for  regions  and  properties  of  interest 
(i.e. ,  those  of  the  data)  the  thermal  properties  of  t ±e  absorbing 
region  are  not  important,  while  those  or  the  surrounding  filar-are.  The 
optical  properties  of  the  absorbing  region,  however,  are  important. 

This  applies  particularly  to  the  imaginary  part  of  die  index  of  refrac¬ 
tion  which  is  a  direct  measure  of  the  energy  that  is  absorbed  in  the 
region. 

Both  observation  {Bliss,  Milam  and  Bradbury  1973]  and  theory 
[Meyer,  Bartoli  and  Kruer  1980]  indicate  that  the  evolution  of  the 
imaginary  part  of  the  index  of  refraction  is  feist  with  respect  to  the 
incident  pulse  length.  Further  theoretical  evidence  of  this  will  be 
presented  in  Chapter  6.  It  is  not  likely  that  these  highly  absorbing 
regions  exist  before  the  incident  pulse. 

In  addition  to  the  effects  of  the  surrounding  film  on  the  laser 
damage  threshold  for  a  single  incident  pulse,  the  mathematical 
generality  of  integral  transform  technique  of  solution  allows  the 
investigation  of  other  effects.  One  such  effect  is  that  of  tne 


substrate  ana  its  conductivity  on  die  rate  at  which  an  absorbing  site 
heats  up.  Another  such  effect  is  that  of  thermal  build  up  due  to  a 
repetitively  pulsed  laser.  The  models  and  results  of  these  studies 
will  be  presented  in  Chapters  3  and  4.  It  is  shown  that  the  substrate 
conductivity  can  be  important  and  that  thermal  build  up  in  repeated 
pulses  is  usually  negligible.  Experimental  coGparison  of  the  work 
described  there  will  be  given  in  Chapter  5.  The  comparison  with 
available  experiments  lends  validity  to  thus  model.  It  also  provides 
an  order  of  magnitude  estimate  of  the  range  of  the  otherwise  unknown 
optical  absorption  coefficient  and  its  behavior. 

Chapter  6  addresses  the  evolution  of  the  electron  density  which 
is  the  source  of  the  imaginary  part  of  the  index  of  refraction-  in  the 
optical  absorption.  It  is  shewn  that  cannon  absorbing  centers  in 
optical  materials  used  for  coatings  can  provide  significant  localized 
absorption  adequate  for  the  evolution  of  the  localized  sites  which  lead 
to  damage.  It  is  also  likely  that  these  centers  exist  in  high  densi¬ 
ties  in  optical  coatings  due  to  the  conditions  of  deposition  and 
polishing. 

Chapter  7  addresses  the  role  of  temperature  dependence  in  the 
absorption  and  diffusion  leading  to  damage  in  optical  films.  The 
findings  are  that  temperature  dependent  absorption  can  (but  doesn't 
necessarily)  contribute  to  damage  in  small  band  asp  oxide  films,  but 
cannot  have  any  significant  effect  in  the  larger  band  gap  fluoride 
films.  Temperature  dependence  of  material  parameters  can  be  important 
quantitatively  for  the  diffusion  processes,  but  the  dependence  is  very 
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caspiicafced  end  often  unknown.  The  dependence,  however,  should  not 
affect  qualitatively  the  findings  of  this  work. 


CHAPTER  2 


THE  LOCALIZED  ABSGRPTIOi-OIFFUSICti  rtOOEL 

The  localized  absorptioo-di f fus i on  model  incorporates  a  snail 
localized  region  within  the  film  that  absorbs  the  incident  laser  radia¬ 
tion  at  a  rate  far  above  that  of  the  surrounding  film.  The  absorbed 
energy  subsequently  diffuses  through  the  surrounding  film.  The  cross 
section  of  the  absorbing  center  available  to  the  radiation  together 
with  the  diffusion  capability  of  the'  film  determine  the  rate  at  which 
the  region  approaches  the  critical  tacperature  at  which  damage  ensues. 

In  this  section  an  analytical  expression  for  the  damage 
threshold  is  described  which  is  valid  for  arbitrary  shapes  of  an 
incident  pulse.  The  particular  example  of  a  single  flat-top-pulse  is 
then  discussed. 

The  derivations  are  shown  in  Appendix  A  and  the  predictions  of 
this  model  are  compared  to  experimental  data  in  Chapter  5. 

A  general  description  of  the  problem  is  illustrated  in  Fig.  1. 

A  spherical  absorbing  region  of  radius  a  is  embedded  in  a  surrounding 
material  of  different  optical  and  thermal  properties.  Only  the  region 
is  permitted  to  absorb  radiation.  The  size  of  the  region  is  assumed  to 
be  ranch  smaller  than  the  laser  spot  size  so  that  absorption  can  be 
thought  of  as  occurring  in  a  field  with  a  constant  radial  intensity. 

The  thermal  properties  of  the  system  are  described  by  the 
thermal  conductivity  K  and  the  thermal  diffusivity  D  of  each  region. 
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Fig.  1.  Quantities  relevant  to  inclus ion-dominated  breakdown  model. 


These  properties  in  concert  with  the  rate  at  which  radiation  is 
absorbed  are  taken  to  be  independent  of  tenperature.  This  does  not 
mean  that  absorption  is  assumed  to  be  constant.  The  assumption  is  that 
the  change  in  absorption  is  sudden  (i.e.,  of  small  duration  cccpared 
with  the  pulse  length) .  The  system  is  said  to  be  damaged  when  the 
interface  betweei  the  two  regions  reaches  seme  critical  temperature 
(e.g. ,  the  melting  temperature  of  the  film).  This  condition  is  chosen 
because  an  interface  between  a  film  and  an  absorbing  region  would  be  a 
stressful  point,  and  melting  has  been  observed. 

Mathematically  the  temperature  distribution  in  the  system  is 
described  by  the  following  thermal  diffusion  equations  [Ozisik  1980] : 
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where  the  subscripts  i  and  h  refer  to  the  inclusion  and  the  host 
respectively  and  A  is  the  source  term.  The  interface  condition  between 
the  impurity  and  the  host  is  assumed  to  be 

Tt  (a,t)  »  Th(a,t)  ,  (4) 

dT  dT. 

tvi  dr  *  ‘Si  dr~  v5) 
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The  solution  for  this  system  may  be  found  vsing  th^  LaPlace 


(temporal)  transform  or  by  an  integral  (spatial)  transform  as  is 
presented  here.  The  integral  transform  technique  is  more  general 
because  no  specific  form  of  the  absorption  profile  need  be  assumed. 

The  general  solution,  for  a  sphere  in  an  infinite  medium,  which 
is  is  derived  in  Appendix  A,  is  found  to  be 


T-(r,t) 


D.b 

l 

2»2raKi 


y2  sin (yr/a)  exp(-y2t/-r)A(y,t) 

2  2  2  2  ’ 
(C  siny  -  y  cosy)  +  b  y  sin  y 


where 


A(y,t)  *  4*  f 


2  l* 

dt'  erpU  t’>  dr  rA{r,t‘)  sin(ar), 
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1/2 

with  y  3  Aa/Di  »  oa.  Now,  if  it  is  assumed  that  A(r,t)  *.-Ar(t) 

(i.e.,  absorption  within  the  region  is  independent  of  the  radius),  and 
A(r,t‘ )  is  proportional  to  the  incident  pulse  which  is  assumed  to  be  a 
square  wave,  a  very  simple  solution  can  be  found. 

Since  the  quantity  of  interest  is  the  energy  of  the  incident 
pulse  which  results  in  an  assuaed  critical  temperature,  the  following 


integrations  are  performed. 


j  dr3A(r,t)  »  Q  ,nJ)I(t) ,  (8) 
vi 

where  q  is  the  absorption  cross  section  computed  from  electromagnetic 
theory  and  n’  is  the  imaginary  part  of  the  index  of  refraction.  The 
dacace  threshold  is  then  defined  as 


f  ?  2 

£  =  j  dtl (t) (J/ca  ) . 

J0 


(9) 


13 


where  t  is  the  pulse  length  required  to  reach  the  critical  temperature 

?c  at  the  radius  a  of  the  absorbing  region.  Thus  Tctp  *»  f{r,Q,E)  is 

inverted  to  give  E  *  g(T  ,t  ,r,Q)  for  the  damage  threshold  in  Joules 

c  p 

per  centixaetar  squared* 

For  a  single  square  pulse  such  that  t^/y  >  1,  where  /  *  a 2/D  , 
in  an  absorbing  region  where  Q  is  proportional  to  Ta2  an  appropriate 
expression  for  the  integral  solution  may  easily  be  computed.  The 
damage  in  Joules  per  centimeter  squared  is  then  found  to  be- 

16T 

e0*^-£(0  hVhV  •  (10) 

This  derivation  is  shown  in  Appendix  A.  Mote  that  the  thermal 
properties  of  the  absorbing  region  do  not  show  up  here.  This  is  an 
approximation  which  is  demonstrated  in  Appendix  A.  Tho  comparison  of 
this  relation  to  experimental  data  is  given  in  Chapter  5. 

It  is  mentioned  above  that  relation  (10)  is  derived  based  upon 
the  fact  that  the  absorption  cross  section  Q  is  proportioned  to  Ta2. 

Although  this  is  usually  thought  of  as  the  large  radius  limit  of  an 
absorbing  sphere,  this  is  not  necessarily  the  case.  A  numerical 
[Kiscaabe  1579]  study  of  the  Hie  theory  of  the  interaction  of  an 
electromagnetic  wave  and  a  generalized  dielectric  sphere  clarifies  this 

i 

point. 

In  general,  the  absorption  cress  section  of  a  sphere  takes  on 
the  appearance  of  Fig.  2.  The  index  of  refraction  of  the  host  material 
and  the  sphere  for  this  case  are  *  1.6,  tu  *  1.9  -  i0.005 
respectively.  The  complicated  structure  results  from  internal  electric 
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and  magnetic  multipole  resonances  which  are  set  up  in  the  interaction 
region.  Energy  considerations  (to  be  discussed  in  Chapter  6)  dictate 
that  for  adequate  absorption  to  cause  damage  within  the  pulse  absorp¬ 
tion  rmist  be  orders  of  magnitude  greater.  Data  fits  of  the  absorption- 
diffusion  model  in  fact  show  that  the  absorption  must  be  such  that  the 
imaginary  part  of  the  index  of  refraction  is  greater  than  one. 

A  study  of  regions  with  an  index  of  refraction  such  as  this 
gives  results  illustrated  in  Fig.  3.  The  important  point  to  note  is 
that  spheres  that  are  only  a  fraction  of  a  wavelength  in  size  tend  to 
behave  as  large  spheres.  This  behavior  results  from  the  fact  that  the 
energy  is  attenuated  within  the  bounds  of  the  sphere. 

As  long  as  the  regions  are  strongly  absorbing  (as  above) 
relation  (10)  bolds.  If  the  region  is  very  stall  compared  to  the  wave¬ 
length,  or  weakly  absorbing,  the  scaling  relation  may  be  shown  to  be 


a 


(ID 


This  relation  has  not  been  shown  to  fit  any  of  the  data. 
Furthermore  the  conditions  under  which  it  holds  true  do  not  allow 
sufficient  energy  to  be  transferred  to  the  region  to  damage  it. 

Relation  (10)  therefore  evolves  frcm  thermodynamic  considerations  of 
energy  transfer  requirements  to  explain  the  observed  morphology  in  some 
of  the  optical  coating  damage  data.  The  morphology  of  pitting  which  is 
the  source  of  this  model  is  primarily  observed  in  the  fluoride  data. 

It  will  be  shown  in  Chapter  5  that  relation  (10)  best  describes  the 
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3.  General  absorption  cross  section  versus  size  parameter  for 
various  imaginary  indices  of  refraction. 


Fig. 


N>|l< 


behavior  of  the  fluoride  data.  The  behavior  of  the  oxide  data  will  be 
explored  more  in  Chapter  7. 

Although  relation  (10)  describes  the  behavior  of  some  of  the 
data  remarkably  well,  it  is  a  proportionality  scaling  of  the  thermal 
properties.  Relation  (10)  says  very  little  about  the  optical  proper¬ 
ties  o tiier  than  the  fact  that  the  imaginary  part  of  the  index  of 
refraction  trust  be  large  for  it  to  be  valid.  Also,  although  m  the 
region  in  which  relation  (10)  is  applied,  it  is  shown  to  be  valid, 
relation  (10'  is  still  an  approximate  solution.  For  these  reasons,  an 
integrated  computer  program  which  solves  the  electromagnetic  field 
equations  for  the  exact  classical  absorption  cross  section  and 
integrates  the  thermal  diffusion  equations  for  the  exact  solutions  of 
them  has  been  developed.  Comparison  of  the  scaling  with  the  exact 
solution  is  shown  in  Figs.  4  through  7.  The  program  is  used  to  find 
the  required  value  that  the  imaginary  part  of  the  index  of  refraction 
oust  assume  in  order  for  damage  to  occur  in  the  films  for  which  damage 
data  exists.  Comparison  of  this  data  with  the  predictions  of  the 
program  is  given  in  Chapter  5.  Hie  scalings  are  shown  to  agree 
remarkably  well  with  the  data.  Hus  point  reinforces  the  concept  that 
very  large  localized  ionization  must  occur  in  a  very  short  time 
compared  with  the  pulse  length  to  explain  the  data.  This  concept  is 
further  reinforced  by  theoretical  computation  carried  out  in  Chapter  6. 
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Table  1.  Quantities  alloyed  in  the  analysis  of  an  absorbing  isoirity 
embedded  in  a  host. 


Thermal  Properties  for  ThF^  (Host)  and  Th02  (Inourity) 
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0.1(J/aa-s-  K) 

K. 
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0.142  (J/crn-s-  K) 

D. 

n 

= 

0.15  (an2/s) 

Di 

- 

0.0594  (an2/s) 

Thermal  Properties  for  Glass  (Host)  and  Platinum  (impurity) 
:<b  =  0.013 (J/cm-s-  K) 

=  0.67  (J/an»s-  K) 

2 

=  0.003  (an  /s) 

Di  -  0.24  (an2s) 

Damage  Condition 
Tc 
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In  Chapter  2  the  dependence  of  the  damage  threshold  on  the 
thermal  and  optical  properties  of  the  system  is  discussed.  The  present 
chapter  specifically  addresses  the  thin  film  aspects.  These  are 
included  by  using  a  cylindrical  impurity  embedded  within  the  film.  A 
description  of  the  model  is  shown  in  Fig.  8. 

The  work  presented  here  is  a  study  of  a  cylindrical  inclusion  ' 
•chadded  in  a  single  layer  thin  film  with  an  mderlying  substrate  whose 
thermal  character  ranges  from  noocooductive  to  highly  conductive.  The 
reason  for  choosing  a  cylindrical  impurity  is  threefold:  cylindrical 
cr  coltnear-llkfi  inclusions  do  occur  in  thin  films  {Lewis,  et  al. 

1985] ,  the  effect  of  the  substrate  is  maximized  (providing  an  upper 
bound)  by  having  a  large  contact  area  (as  opposed  to  a  point  contact 
such  as  with  a  sphere) ;  and  finally,  it  affords  a  mathematical 
simplicity  through  the  latching  of  boundary  conditions. 

In  reference  to  Fig.  8  a  cylindrical  impurity  of  radius  a  is 
embedded  in  a  thin  film  of  physical  thickness  1  with  different  thermal 
and  optical  properties.  The  diameter  of  the  inclusion  is  assumed  to  be 
much  less  than  the  spot  size  of  the  laser  illuminating  the  sample  sc 
that  the  impurity  is  subject  to  a  uniform  intensity  in  the  radial 
direction.  The  film  and  absorbing  region  are  in  perfect  thermal 

contact  with  an  underlying  substrate.  The  thermal  conductivity 

-1  -1  - 1.  2  -1 
K(J  aa  s  K  ;  and  the  thermal  aiffusivity  D(cm  -s  )  are  labeled  by 
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Pig.  8.  A  Description  of  the  General  Film-Substrate  Model 

A  cylindrical  impurity  of  radius  a  is  embedded  in  a  thin  film 
of  thickness  1  ana  interfaces  with  a  substrate.  The  thermal 
properties  of  tbs  impurity/  host  and  substrate  are  subscripted 
with  i,  h,  and  s  respectively.  The  diameter  of  the  inclusion 
is  -assaaed  to  be  much  scalier  than  the  scot  size  of  the  laser 
illuminating  the  sample.  Radiation  is  absorbed  in  the 
inclusion  and  the  diffusion  governs  the  tsnperature  field  that 
results. 
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the  subscript  i,  h,  or  s  referring  to  the  inpurity,  film  (host)  and 
substrata  respectively.  The  thermal  properties  of  the  system  along 
with  the  absorption  cross  section  are  taken  to  be  independent  of 
temperature  and  the  impurity  is  the  only  part  of  the  system  that 
absorbs  incident  radiation.  As  is  the  case  for  the  previously  studied 
spherical  model,  damage  is  assured  to  occur  when  the  surface  of  the 
inpurity  reaches  same  critical  material-dependent  temperature  Tc« 

Mathematically  the  model  is  described  by  the  following  thermal 
diffusion  equations  [Ozisik  1980} : 
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where  A  is  the  source  term  in  watts  per  cubic  centimeter. 
The  boundary  conditions  are  assuned  to  be 
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d?h  d? 
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(18) 


and 

||  *  0  at  2  -  Q  .  (19) 

Exact  solutions  are  found  for  the  extreme  cases  of  a 
nonconducting  substrate  and  an  infinitely  conducting  substrate  with 
either  uniform  or  standing  wave  absorption  distributions.  As  before 
the  solution  may  be  found  via  Laplace  or  integral  transform.  An 
approximate  solution  is  also  found  for  the  important  realistic  case  of 
a  finite  valued  thermally  conducting  substrate.  All  of  these  cases  are 
compared  with  cuaerical  finite  element  analyses  via  the  thermal/ 
structural  analysis  code  KASTRAU  for  consistency.  The  derivation  of 
the  results  are  presented  in  Appendix  B.  The  Laplace  derivation  is 
raenCicoed,  but  not  presented  due  to  its  extraordinary  length  and  its 
lade  of  generality.  It  does  not  apply  to  arbitrary  absorption 
functions  in  space  os  time.  It  is  found  oscsiAtcet  with  t h»  integral 
transform  solution. 

For  the  special  cases  of  the  infinite  and  nonconducting 
substrates  the  problem  reduces  to  Eqs.  (12)  and  (13)  together  with  Bq. 
(15)  '-her*  Eqs.  (17)  and  (18)  reduce  toT*0atz*l  for  an 
infinitely  conducting  substrata  and  dT/as  *  0  at  z  »  1  for  a 
nonconducting  substrate. 

The  general  solution  for  the  boundary  temperature  is  presented 
in  Appendix  3.  Any  arbitrary  spatial  dependence  for  absorption  may  be 
used.  Any  such  function  tray  also  be  repetitively  pulsed  as  is  shown  in 
the  next  chapter. 


For  the  purpose  of  demonstrating  trends/  the  realistic  case  of 


a  single-pulsed  sine  squared  standing  wave  (half-wave  film)  is 
investigated.  The  source  term  say  be  written  as 
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The  source  term  A(zrt>  can  also  be  written  in  terns  of  the 
incident  laser  intensity  ICKcas ”'i)  and  the  absorption  cross  section 
C(an2)  by  realizing  that 
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where  tp  is 


the  duration  of  the  laser  pulse. 
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In  this  case 
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As  before,  the  above  relation  is  integrated  ovts:  t  and  inverted 


to  ?it»d 


E(J/ca2)  »  g(Ti?  tp,  r,Qj .  (24) 

The  absorption  cross  section  used  in  this  solution  is  u  *  Cic~ 
where  C  is  sea*  constant.  This  is  based  upon  the  Mie  scattering 
results  from  Chapter  2. 

Once  the  solutions  for  the  damage  threshold  (in  Joules  per 
square  centimeter)  for  a  given  pulse  length  are  known,  the  thin  film 
character  and  substrate  effects  in  extreme  limiting  cases  can  be 
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determined.  In  order  to  isolate  the  substrate  effects  from  radial 
diffusion,  the  limit  as  a  »  is  taken.  Expressions  for  the  damage 
threshold  scaling  are  found,  for  the  case  of  K$  -*•  0,  to  be 


2D  ;t  •»  -  2  1 

E  5  V^i  (~i —  +  -^2  H  -  exp  l  -  {—)  Dit]})'1  , 
4~  i 


(25) 


and  furthermore  as  t  -  0 
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For  the  long  pulse  limit,  i.e.  t  s>  1/(4D^~  )  we  find 
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3y  Caspar ison ,  for  a  substrate  which  is  infinitely  conducting. 


as  t  -  0 
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T  IK. 
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which  is  the  same  as  is  found  for  the  nonconducting  substrate  above. 
Finally,  for  the  long  pulse  lioit  with  an  infinitely  conducting 
substrate,  i.e.  t  »  412/(9&i*2) , 
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The  important  scaling  trends  of  the  damage  threshold  are  the 
linear  dependence  on  temperature  T  ,  the  linear  dependence  on  the 
thickness  1  of  the  film  for  very  short  pulse  lengths  (t  ■*  0) ,  the 
inverse  thickness  dependence  on  1  for  long  pulses  with  an  infinitely 
conducting  substrate,  and  the  various  pulse  length  dependences  (e.g. 
none  as  fc  -  0  and  linear  in  t  for  long  pulses  with  infinitely  conduct¬ 
ing  substrates) .  Some  of  these  dependences  begin  to  appear  in  Fig.  9 
for  the  range  of  pulse  lengths  plotted.  For  example,  we  begin  to  see 
an  inverse  thickness  dependence  for  long  pulses  and  a  highly  conducting 
substrate. 

It  should  be  noted  that  these  scalings  are  in  extreme  limits 
and  are  not  intended  to  be  realistic  design  relations.  For  example,  no 
host  dependence  appears  because  of  the  limit  a  -  These  relations 
are  intended  to  give  insight  into  the  effects  or  lack  of  effects  of 
thin  film/ substrate  system  parameters. 

The  general  thermal  parameter  dependences  of  the  problem  due  to 
radial  heat  diffusion  are  determined  by  insulating  the  substrate  and 
varying  the  host  film  and  inpurity  thermal  conductivities.  As  with 
previous  spherical  impurity  studies,  the  damage  threshold  increases 
with  increasing  host  thermal  conductivity,  while  the  carnage  threshold 
is  essentially  insensitive  to  the  thermal  conductivity  of  the 
inclusion. 

A  more  illuninatiog  exercise  is  to  set  1  *  (4/3) a  so  that  the 
length  of  the  inclusion  varies  with  the  radius  of  the  inclusion.  This 
forces  the  volisne  to  be  proportional  to  'a*  cubed  and  the 
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cross  section  to  be  proportional  to  ’a’  squared.  When  this  is  done  a 
minimum  value  of  the  damage  threshold  is  found  {i.e.  a  radius  at  which 
damage  is  easiest)  and  scaling  laws  for  host  conductivity  and  pulse 
length  are  found  much  like  those  predicted  for  spherical  impurity 
particles.  The  end  result  of  this  analysis  is  that  the  details  of  the 
geometry  are  not  very  important  to  the  thermal  aspects  of  the  problem 
as  long  as  the  volune  of  the  impurity  is  proportional  to  seme  mean 
radius  cubed  and  the  absorption  cross  section  is  proportional  to  some 
mean  radius  squared.  So,  if  the  absorbing  regions  are  irregularly 
shaped  volumes,  as  may  be  the  case,  then  the  thermal  scaling  laws 
derived  from  an  ideal  geometry  model  should  still  apply. 

The  final  aspect  of  the  problem  studied  is  the  effect. -of  the 
intensity  distribution  within  the  film.  It  is  easily  seen  from  Fig.  10 
that  high  intensity  I  (J  cm~2s~i)  short  pulses  can  lead  to  very  high 
temperature  gradients.  If  the  intensity  is  high  (Fig.  10,  curve  3) 
with  a  sine  squared  half-wave,  and  the  pulse  length  short  enough,  very 


high  temperatures  can  be  reached  in  the  center  of  the  film  while  the 
boundaries  remain  nearly  at  their  initial  temperature.  In  this  case  a 
substrate  will  have  very  little  effect  regardless  of  how  highly 
conducting  it  may  be. 

The  descriptors  short,  long,  etc.  refer  to  a  comparison  with 
scales  determined  by  the  diffusivities  D  involved.  For  example,  the 
diffusion  length  5  »  (D  t)  for  I.  is  approximately  1/40  of  the 
thickness  of  the  film.  At  this  depth  in  the  film  the  temperature  is 
still  quite  low.  Thus,  this  rough  measure  of  the  range  of  the  effect 
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of  one  point  or.  another  indicates  that  the  boundary  temperature  will 
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Plot  of  Temperature  vs.  Film  Depth  -  I. 

The  absorbed  energy  is  proportional  to  a  standing  wave  for 
three  different  intensities  Ilf  I2,  I3  (I3  >  l2  >  I 1 ) ,  three 
pulse  lengths  and  a  fixed  energy  density  Z  =  0.9  J  an'-: 
curve  1,  1 1 1  tp  =  10-7  s;  curv’e  2,  I,,  tp  =  10~3  s;  curve  3, 
1 1,  tp  =  10" '  s.  Tbs  sane  total  energy  is  deposited  m  all 
three  cases  (i.e.  It  =  constant). 
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not  be  affected  significantly.  However,  if  we  decrease  the  thickness 
of  the  film  by  a  factor  of  ten  with  the  same  pulse  length,  6  is 
approximately  1/4  of  the  film  thickness.  Referring  to  Fig.  11  we  see  a 
significant  effect  or.  the  temperature  profile  of  substrates  with 
various  conductivities. 

If  instead  of  a  half-wave  intensity  distribution,  we  have  a 
wave  such  that  the  saxirrarn  intensity  is  at  the  interface  with  the  film 
and  substrate,  there  would  be  a  considerable  effect  for  good  conducting 
substrates  at  much  smaller  <5  parameters. 

The  results  presented  here  imply  various  conclusions.  At  least 
one  of  the  conclusions  presented  is  easily  verified.  That  is,  if 
substrates  of  high  thermal  conductivity  are  used,  the  damage  threshold 
nay  increase.  Experimental  results  pertinent  to  this  are  presented  in 
Chapter  5. 
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Fig.  11.  Plot  of  Temperature  vs.  Film  Depth  -  II. 

The  absorbed  energy  is  proportional  to  a  standing  wave 
showing  the  effect  of  three  different  substrate  thermal 
conductivities  (t  *  10' -  s,  half-wave  standing  wave 
absorption  profile,  fixed  energy  density  (0.09  J  enf 2) 
TfcF  host;  T11O2  inclusion) :  curve  1,  Ks  =  0  J  cn'^s'1 
(free-standing  film);  curve  2,  Ks  =  0.05  J  cn'!s'lK"’1 
(fused  silica);  carve  3,  Kj  =  3.7  J  cm' *  s' ’•  K~ ‘  (copper) 
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CHAPTER  4 


REPETITIVELY  PULSED  LASER  COtSIDERATICiS 

The  problem  of  multiple  poise  laser  induced  damage  in  optical 
components  is  a  very  inportant  and  general  one.  It  is  a  significant 
problem  that  reflects  the  durability  of  optical  elements.  This  stars 
from  the  fact  that  multiple  pulsed  dasage  addresses  the  use  of  optical 
elements  that  are  subjected  to  more  than  one  cycle  regardless  of  the 
time  delay  between  pulses. 

In  a  review  of  this  area  one  encounters  a  large  array  of 
seemingly  contradictory  data.  For  example  in  same  instances  -there  is 
no  apparent  effect  of  prior  pulses  [Bass  and  Barrett  1972]  (i.e.,  no 
memory) ,  while  in  other  cases  there  is  unmistakable  evidence  of  a 
build-up  or  accunulation  of  damage  [Merkle,  Bass  and  Swim a  1982j .  Seme 
situations  demonstrate  clear  evidence  of  thermal  build-up  effects 
[Kounvakalis,  Lee  and  Bass  1982;  Wood,  Sharma  and  Wate  1982]  while 
others  [Bordelon,  et  al.]  show  no  evidence  of  this  effect.  A  pulse 
repetition  frequency  dependence  in  scrae  cases  [Kouavakalis,  et  al. 

1982;  Wood,  et  al.  1982;  Balitskas  and  Maldutis  1981]  exhibits  revers¬ 
ible  effects  in  some  materials  while  others  [Eordeloo,  et  al  1982]  show 
no  evidence  of  this  effect.  It  is  also  known  that  low  level  irradia¬ 
tion  can  condition  elements  such  that  they  can  ultimately  withstand 
higher  fineness  [Gorshkov,  et  al  1983;  Becker,  et  al.  1983).  This  wide 
diversity  of  responses  may  also  depend  upon  the  configuration  of  the 
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material  (i.e.  bulk-  film,  crystalline,  amorphous,  aetallic,  polymer, 
etc.) 


The  morphology  of  nail ti -pulsed  damage  is  just  as  varied,  not 
only  between  different  materials,  but  between  different  configurations 
of  tbe  same  material  (i.e.,  anti-reflecting  films  and  highly-reflecting 
films)  (Foityn  and  Jolin  1983] .  One  mode  of  damage  foi  films,  both 
with  and  without  setae  sort  of  accumulation  mechanism,  is  for  a  small 
pit  to  form  and  no  subsequent  damage  to  occur  [Fcltyn  and  Jolin  1983] . 
Another  mode  is  for  damage  to  grow  with  subsequent  pulses  to  a  diameter 
that  is  consistent  with  points  in  the  incident  field  that  are  *v  25%  of 
tbe  damage  threshold  [Fcltyn  and  Jolin  1983] .  Clearly  there  is  no 
sisple  universal  multi  pulse  laser  induced  damage  scenario.  There,  in 
fact,  may  be  as  many  effects  as  there  are  materials  or  classes  of 
materials. 

What  is  done  herein  is  to  once  again  exploit  the  thermal  model 
of  laser  induced  damage  and  apply  it  to  multi  pulsed  lasers.  The  justi¬ 
fication  for  this  is  that  absorption  and  diffusion  must  occur  in  all 
cases  of  laser  damage  and  laser  induced  damage  is  thermally  dominated 
in  certain  single  pulse  cases.  There  is,  in  addition,  evidence  that 
the  mechanism  for  some  cases  of  laser  induced  damage  is  the  same  for 
tbe  first  pulse  as  the  subsequent  pulses  [Balitskas  and  Kaldutis  1981; 
Bass  and  Barrett  1973;  Kovalev,  at  al.  1980]  (though,  for  other  cases 
(Kounvakaiis,  et  al,  1982]  this  is  evidently  rot  true).  Even  cases 
where  damage  is  not  thermally  dominated  a  change  of  state  necessarily 
accompanies  a  change  in  temperature  (e.g.  a  change  in  band  gap)  and 
thus,  it  is  desirable  to  knew  tbe  temperature  field. 
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The  thermal  model  is  generalized  to  facilitate  any  sort  of 
general  absorption  that  may  occur  in  various  materials.  However,  the 
absorption  function  must  be  determined  for  each  different  situation. 

The  determination  of  fee  temporal  and  spatial  dependefjce  of  fee 
absorption  function  for  each  situation,  as  for  the  single  pulse  case, 
is  fee  primary  difficulty.  For  certain  single  pulse  situations  a  very 
simple  absorption  function  can  prove  adequate.  However,  due  to 
irreversible  changes  (i.e.  manory}  of  seme  sort,  this  is  not  true  for 
most  multiple  pulse  interactions. 

In  this  section  a  general  model  of  arbitrary  multi-pulse 
absorption  is  derived.  The  arbitrary  form  of  the  absorption  applies 
not  only  to  localized  a^orption,  but  allows  application  to  laser 
windows  with  macroscopic  gaussian  beam  absorption  [Mood,  et  al.  1982} « 
In  addition,  it  allows  application  to  an  opaque  material  with  absorp¬ 
tion  of  any  general  form.  'She  problem  is  then  simplified  to  the  case 
of  a  uniform,  nonexpanding  absorption,  without  irreversible  processes 
as  an  illustration  of  fee  effects  purely  due  to  tberrral  build-up. 

The  results  are  as  expected,  in  that  thermal  build-up  in  itself 
cannot  be  responsible  for  the  observed  decrease  in  damage  threshold  in 
most  cases  of  isolated  absorbing  inclusions.  Other  mechanises  must 
couple  with  the  diffusion  process  to  provide  the  observed  behavior. 

The  key  to  understanding  fee  many  observations  and  mechanisms  involved 
is  to  determine  specific  pulse  to  pulse  behavior  of  fee  absorption 
function.  Any  correct  theory  of  fee  interaction  involved  must  predict 
the  spatial,  temporal,  pulse  repetition  frequency  and  pulse  number 
dependence  of  the  absorption  term  in  the  diffusion  equations. 
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Theory 

In  previous  chapters  a  spherical  model,  valid  for  the  region 
where  the  pulse  length  is  not  larger  than  the  diffusion  time  to  the 
boundaries  of  the  film,  and  a  cylindrical  model,  valid  for  ouch  longer 
times,  are  described. 

In  general  the  diffusion  equation  is  used  to  describe  the 
distribution  of  temperature  within  the  system.  In  the  absorbing  region 
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where  T,  D,  K  are  temperature,  diffusivity  (ca2/s)  and  conductivity 
(J/  K-s-cn)  and  A  is  the  absorption  fraction  in  (J/cm^-s)  while 
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outside  of  the  absorbing  region.  Here,  perfect  thermal  contact  is 
assumed  at  the  boundaries  of  the  absorbing  region  and  the  subscripts  i 
and  h  refer  to  the  inclusion  and  host  film  respectively.  The  two 
geometries  give  specific  solutions  within  the  impurity  of 
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Hers  A  is  the  transformed  absorption  function  (defined  below) , 
N(a)  is  a  normalization  factor,  a  the  radial  eigenvalue  and  1  is  the 
length  of  the  cylinder  (i.e.  the  thickness  of  the  film).  The 
expression  for  the  source  term  in  the  above  equations  is 
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Here  5  is  the  composite  eigenvalue  of  the  system,  (r,a)  is  the 
eigenfunction  of  the  composite  region  and  A(r,t)  is  tee  general 
absorption  function  with  *Vola  referring  to  tee  volume  of  the  absorbing 
region.  The  subscript  j  on  t>  applies  to  axial  modes  of  the  cylindrical 
geometry.  Expressions  for  these  parameters  are  given  in  Table  2  for 
specific  geometries.  The  absorption  function  written  in  the  form 
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contains  the  complete  evolution  of  the  absorption  including 
irreversible  or  long  lived  changes  to  the  material  with  each  pulse  from 
1  to  N.  There  is  data  (Bliss,  et  ai  1973;  Bordelon-  et  al  1982) 
indicating  that  in  some  cases  tee  change  in  absorption  occurs  within  a 
time  teat  is  small  compared  to  tee  delay  time  *td*  between  pulses.  In 
some  cases  it  may  occur  within  a  time  that  is  snail  compared  to  the 
pulse  length  itself.  For  an  idealized  case  where  the  absorption 
evolves  early  into  the  pulse  and  tee  absorption  teat  follows  is 
proportional  to  the  field,  or  the  case  where  the  evolution  transient 
occurs  completely  in  between  the  puises 
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Table  2.  Eigenfunction  parameters  for  a  spherical  and  a  cylindrical 
inclusion 
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A(r,t)  *  l  ^Wlsit-ln-Dt 

* 

*  (n-l)td3)-e((t-intp  v  (n-Dtj])}.  (26) 

Here  the  nth  pulse  tarns  on  at  t  3  {n-l)t  +  (n-l)t,  and  turns  off  at 

p  d 

t  3  ntp  (R-l)t^  with  t^  being  the  pulse  length  and  tj  being  the  delay 
time  between  pulses  and  3  being  the  step  function. 

ttnxs,  the  spatial  dependence  of  absorption  from  pulse  to  pulse 
is  all  that  is  important.  In  other  cases  [Foltyn  and  Jolin  19331 
nul ti pulse  damage  has  been  observed  to  occur  after  some  number  of 
pulses,  but  not  grow  spatially,  in  fact,  no  further  damage  is  observed 
after  the  initial  pitting.  In  this  case 

&(r,t)  »  f (r)  l  A  t 6 (t- [ (n-1) c 

DU  P 

+  (a-DtdJ)  -  9(t-[ntp+  (n-l)tdl)}  ,  (37) 

so  that  only  the  magnitude  of  absorption  varies  with  the  pulse  number. 

In  all  of  these  cases  the  absorption  evolution  is  unknown  and 
is  observed  to  be  material  and  laser  parameter  dependent  with  an 
unknown  physical  mechanism. 

'Whatever  the  driving  source,  a  measured  empirical  or 
theoretically  derived  fora  of  A  most  be  found  before  any  specific 
situation  can  be  treated  adequately. 

An  illuai  rating  exercise  as  a  first  step,  however,  is  to  assume 
A  is  constant  and  develop  the  solution.  There  are  sane  cases  where 
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Table  3.  Transform  coefficients  for  various  absorbing  regions. 
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pulsed  laser  induced  damage,  it  is  likely  that  it  does  apply  to  some 
cases  and  may  be  informative  for  others  for  which  it  is  not  strictly 
valid.  It  is  informative  because  it  illustrates  the  purely  thermal 
effects. 

Results 

Nucbers  are  applied  to  the  previous  expression  providing  a 
sense  of  the  effect  of  thermal  build-up  in  optical  materials  subjected 
to  repetitive  laser  pulses.  For  exaiple,  for  a  one-half  micron  thick 
thoriixn  fluoride  film  using  the  isolated  spherical  absorbing  impurity 
model  with  a  0.1  nanosecond  pulse,  if  no  growth  of  the  absorbing  region 

c 

occurs  a  repetition  frequency  of  greater  than  10  Hz  would  be  needed  to 
see  any  significant  decrease  in  damage  threshold  due  to  thermal  build¬ 
up.  live  criterion  that  is  used  for  damage  there  is  for  the  edge  of  the 
absorbing  region  to  reach  the  film  melting  temperature.  Note  on  Fig. 

12  the  increase  in  radius  most  likely  to  damage  with  increased  number 
of  pulses.  So,  if  no  growth  of  the  inclusion  occurs  (i.e.  disc  like 
growth  into  the  film)  and  growth  is  limited  by  the  film  thickness 
then  the  effect  of  thermal  build-up  is  even  less.  Note  also  in  Fig.  13 

C 

of  pulse  repetition  frequency  equal  to  10  Hz,  that  essentially  no 
thermal  build-up  occurs.  Thus  thermal  build-up  could  not  be  a  factor 
in  damage  for  this  idealized  case.  This  does  not  take  into  account  a 
possibly  dense  distribution  of  particles  so  that  the  diffused  heat  of 
nearby  absorbing  regions  could  build-up  in  the  surrounded  region. 
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Fig.  12.  Coeluted  damage  threshold  in  {3 /an2)  versus  radius  of  the 
absorbing  inclusion  -  I 

This  is  a  plot  of  one  example  of  computed  damage  threshold  in 
(J/an2)  versus  radius  of  the  absorbing  inclusion  for  various 
numbers  of  pulses.  The  pulse  repetition  frequency  is  10-  Hz. 
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13. 


Ccccutad  damage  threshold  in  (J/crs2)  versus  radius  of  the 
absorbing  inclusion  -  II 


Ibis  is  a  plot  of  caroutad  damage  threshold  in  {J/car)  versus 
radius  of  the  absorbing  inclusion  similar  to  Fig.  1,  but  vitn 
a  pulse  repetition  frequency  of  105  Hz. 


Detailed  experimental  determination  of  the  distribution  and  strength  of 
absorbers  would  be  required  for  this  calculation. 

The  largest  possible  effect  of  thermal  build-up  for  a  single 
isolated  absorber  would  be  for  a  free  standing  film  (allowing  only  two 
dimensional  neat  diffusion)  of  very  lew  thermal  conductivity  with  a 
large  strongly  absorbing  inclusion  and  a  long  pulse  length.  So,  asstxne 
a  free  standing  film  of  Si02  (k  =  0.014  J/c a-  K-s)  tea  microns  thick 
with  a  10  microsecond  incident  pulse  repeated  at  1000  Hz.  The  largest 
inclusion  that  could  initially  be  embedded  in  the  film  would  be  one 
with  a  diameter  equal  to  the  film  thickness,  again,  assuning  no  growth 
of  the  absorbing  region  a  decrease  in  damage  threshold  of  only  ten 
percent  for  1000  pulses  due  to  thermal  build-up  occurs.  For  a- 
substrate  of  greater  than  zero  conductivity  this  value  would  be  further 
decreased.  These  extreme  conditions  should  indicate  that  thermal 
build-up  by  itself  is  not  a  factor  in  these  models.  However,  there  are 
conditions  with  and  without  the  assistance  of  additional  mechanisms 
when  it  is  imaginable  that  thermal  build-up  would  be  a  factor. 

Suffice  it  to  say,  a  dense  distribution  of  absorbing 
impurities,  irreversible  changes  or  ceoperature  dependent  absorption 
that  is  initiated  at  an  absorbing  inclusion  are  a  few  such  cases.  The 
latter  suggestion  is  not  a  new  idea.  Kamolov  [1982}  has  investigated 
noolinaar  temperature  dependent  absorption  cf  the  form 

X(T)  =*  ti  ♦  .  (41) 

Where  X  is  the  impurity  absorption  and  X^e  U/T  is  an  intrinsic  film 
absorption.  This  is  addressed  in  Chapter  7. 


Also  M2nenkcv,  et  al.  [?9S2]  have  suggested  irreversible 


mechano-chemical  resctiore?  with  a  rate  constant  of 
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(44) 


where  u  is  an  activation  energy  and  is  the  stress  proportional 
factor  driven  by  temperature  gradients  assisting  in  the  damage. 

With  one  of  these  mechanists  the  region  that  is-  absorbing 
expards  into  the  film  from  its  origin  at  an  inclusion  or  impurity. 
Another  such  process  is  investigated  in  Chapter  6.  This  expansion  to 
macroscopic  dimensions  then  necessarily  involves  thermal  build-up.  In 
fact  according  to  the  theories  mentioned  above,  the  expansion  in  turn 
is  driven  by  the  thermal  build-up. 

Another  situation  involving  thermal  build-up  is  the  case  where 
absorption  occurs  on  a  macroscopic  scale  (on  the  order  of  the  diffusion 
length  between  pulses) .  This  occurs,  for  example,  in  semiconductor 
grade  windows  [Wood,  et  al  1982]  and  on  metallic  mirrors  [Koumaakalis, 
et  al.  1982].  Table  2  lists  the  spatial  part  of  the  solutions  for  scene 
of  these  cases.  In  particular  the  case  of  a  weakly  absorbing  window 
such  that  absorption  is  uniform  with  depth  into  the  window  and  the  case 
of  a  standing  wave  in  a  thin  film  are  shown.  Many  other  similar 
solutions  are  possible  (surface  absorption  or  absorption  which  has 
essentially  any  arbitrary  depth  dependence) .  These  are  just  mentioned 
to  point  out  the  trenendous  generality  of  the  integral  transform 
technique. 
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As  with  microscopic  absorption,  most  cases  of  observed 
macroscopic  absorption  observed  are  also  dependent  on  tbs  mirbsr  of 
pulses  and/or  the  pulse  repetition  frequency.  So  once  again  a  function 
versus  time,  temperature,  pulse  number  and  pulse  repetition  frequency 
must  be  known  or  assured  to  properly  model  this  situation. 

Some  work  by  Nathan,  felkar  and  Guenther  [1983]  has  indicated 
that  a  multi-pulse  oarage  decrease  can  occur  by  thermal  build-up  and 
various  relations  involving  multipulse  parameters  are  realized. 

Together  these  lead  to 
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(45) 


The  terms  above  are  defined  as  follows:  Ej  is  the  damage 

threshold  in  (J/cm2)  and  N  -  is  the  number  of  pulses  at  which  damage 

sin 

is  aost  likely  to  occur  given  some  pulse  length  and  frequency.  As  in 

Chapter  2  this  was  found  baser*  upon  the  fact  that  the  heat  balance 

gives  seme  particle  radius  rrost  likely  to  damage.  Here  is 

the  radius  of  the  particle  most  likely  to  damage  at  this  pulse  length 

and  pulse  repetition  frequency.  These  authors  mention  that  relation 

(45)  also  applies  to  uniformly  absorbing  semi-infinite  surfaces. 

The  naaerical  work  of  Nathan,  et  al.  [1983]  was  verified  and 

the  basic  trends  are  shown  in  Fig.  12.  In  addition,  approximate 

relations  for  E  (K,o ,C  ,t  ,t .,N)  and  a  .  (K,p ,C  .t  ,t,,N)  *«re  derived 

v  p  c  axn  v  p  a 

exactly  as  in  Chapter  2  further  verifying  the  miaerics.  Unfortunately, 


as  also  indicated  in  Fig.  12,  the  miniisxa  radius  for  damage  to  occur 
quickly  becomes  very  large.  In  fact,  it  becomes  several  orders  of 


magnitude  larger  than  most  thin  films  for  reasonable  laser  parameters. 
Nathan's  relations  therefore  primarily  apply  to  bulk  materials  with 
extremely  large  inclusions,  or  more  reasonably  to  surface  absorption  on 
a  seed- infinite  absorber. 

The  derived  relations  mentioned  above  are  not  given  because  of 
their  extreae  ccnplexity  yielding  little  additional  insight  into  the 
problem.  However,  as  with  most  real  lasers,  if  one  assumes  that  td  » 
t^and  that  impurity  radii  are  snail  (limited  to  a  thin  film  thickness) 
the  relations  simplify  to  the  ones  given  in  Chapter  2  for  the  single 
pulse  case.  This  simply  indicates  that  no  thermal  build-up  can  occur 
(i.e.  in  thin  films,  if  inclusions  damage  thermally,  it  must  be  on  the 
first  shot  instead  of  the  Nth) .  It  also  indicates  that  multi  poise 
thermal  damage  scales  against  the  thermal  parameters  in  a  similar 
manner  to  the  single  pulse  case.  This  has  noc  beers  investigated  for 
pulse  dependent  absorption. 

The  results  of  this  chapter  are  certainly  not  conclusive.  The 
problem  is  as  initially  stated  that  there  is  no  simple  universal  effect 
in  multi puls*  damage  that  we  have  ascertained. 


CHAPTER  5 


COMPARISON  WITH  EXPERIMENTAL  DATA 

the  ranker  of  quality  sets  of  experimental  laser  induced  damage 
of  optical  coatings  data  are  few.  Comparison  of  data  points  frctn 
different  data  sets  is  probably  not  valid.  This  is  due  to  the  range  of 
incomparable  techniques,  differences  in  conditions  and  criteria  for 
damage  and  differences  in  coating  sets.  Thus,  it  is  necessary  to  have 
large  sets  of  parameter  variation  in  a  given  experi-  meet.  There  is 
one  particularly  thorough  set  of  data  in  the  literature  at  present 
[Waiter,  et  al  1981]. 

This  data  set  is  compared  with  the  scaling  relation  derived  in 
Chapter  2  for  similar  films.  The  term  similar  films  refers  to  the 
carnations  of  a  half-wave  film  with  a  half-wave  film,  a  full-wave  film 
with  a  full-wave  film,  etc.  The  reason  for  this  similarity  is  that  the 
wavelength  dependence  in  Mie  theory  appears  in  the  form  of  a  size 
parameter  q  *  2-a/X.  Where  a  is  the  radius  of  the  absorbing  region. 
Although  the  radius  of  the  absorbing  region  is  of  a  size  most  likely  to 
be  danaged  according  to  thermodynamic  considerations,  for  very  thin 
films  it  is  roughly  determined  by  the  thickness  of  the  film.  Thus,  if 
the  film  thickness  is  halved  and  tbs  incident  wavelength  is  halved  the 
size  parameter  q  regains  the  same  and  so  should  the  ratio  of  the 
absorption  cross  section  to  the  geometric  cress  section.  This  in  fact 
is  not  observed.  A  wavelength  dependence  still  retains  in  the 
experimental  data. 


When  the  data  are  plotted  against  the  scaling  (Eq.  10)  for 
half-wave ,  quarter-wave  and  full-wave  films,  a  clear  splitting  by 


wavelength  occurs  in  each  of  the  three  types  of  film.  The  wavelength 
dependence  can  be  realized  by  comparing  Figs.  14-16  with  Figs.  17-19. 
The  only  free  parameter  between  Figs.  17-19  is  the  imaginary  part  n*  of 
the  index  or  refraction.  That  is,  a  value  of  n*  is  chosen  that 
normalizes  the  analytical  curves  to  an  experimental  data  point.  Since 
the  absorption  (1/cmj  *  4sn'/X,  the  nonlinear  dependence  of  n'  on 
indicates  that  the  absorption  is  a  function  of  the  wavelength.  This  is 
not  unreasonable.  Absorption  is  a  function  of  wavelength  in  most 
theories  of  absorption  that  may  apply  here.  The  specifics  of  the 
wavelength  dependence  cannot  be  hypothesized  though  unless  precise 
mechanism  or  mechaniaas  are  known.  It  should  be  noted  that  the 
imaginary  part  of  the  index  of  refraction  is  simply  a  phenomenological 
measure  of  absorption  and  does  not  give  any  information  cn  the 
absorption  mechanism. 

For  oxides  (Fig.  20)  there  is  a  splitting  of  the  scaling  for 
both  wavelength  and  material,  although  HfO-,  scales  on  the  same  line  as 
Al.,0-.  If  oxides  are  scaled  against  t  ^^(KcC^  ^ L  multiplied  by 

i /*> 

same  additional  material-dependent  term  (other  than  (KcCv)  the 
oxides  also  scale  as  in  Eq.  (10).  According  to  data  (Walker,  et  al. 
1981)  this  constant  is  highly  material  dependent.  It  does  not  seen  to 
depend  upon  the  film  thickness,  but  it  does  appear  to  depend  on  the 
incident  wavelength  to  seme  fractional  newer  (about  1/2),  i.e.  the 
magnitude  of  the  splitting  of  each  material  at  each  wavelength  varies 
roughly  as  Ail  this  indicates  that  oxides  also  scale  as  the 
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Scaling  of  the  «xp*x imen  tally  measured  damage  threshold  of 
fluorides  vs.  a  theoretically  derived  parameter  -  I. 

The  garaoeter  contains  the  thermal  properties  and  the  data 
is  for  qaartas-wave  films  at  two  different  wavelengths  (the 
anomalously  high  points  should  be  noted)  measured  at 

5  ns  (I)  and  15  ns  ( 1  ):  -  ,  \  *  1.06  microns; - , 

\  *  0.53  microns.  The  material  properties  for  the  work  m 
this  chapter  come  from  Goldsmith  [1361],  Sparks  [1977]  and 
Touioukain  (1970].  Room  Temperature  values  are  chosen. 
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Fig.  IS.  Scaling  of  the  experimentally  measured  daeage  threshold  of 
fluorides  vs.  a  theoretically  derived  parameter  -  II. 

The  parameter  contains  the  thermal  properties  and  the  data  is 
for  half-wave  files  at  four  different  wavelengths  measured  at 

5  ns  (!)  and  15  ns  (  f  ):  -  ,  \  *  1.06  microns; - , 

x  =  0.53  microns; - ,  \  -  0.35  microns;  . ,  ;  =  0.26 

microns. 
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Fig.  15,  Scaling  cf  the  exp*ximen tally  treasured  dec sage  threshold  of 
fluorides  vs.  a  theoretically  oerived  paran^t-or  -  III. 

The  par  acre  ter  contains  the  thermal  properties  and  the  data  is 
for  full-wave  files  at  three  different  taveiengths  measured 

at  5  ns  (I)  and  IS  ns  (  *):  - ,  >  =  1.06  microns; - , 

X  =  0.53  microns;  . ,  \  =  0.26  microns. 


FlUt  THICKNESS  (X  ) 


Damage  threshold  data  and  theory  vs.  film  thickness  -  II. 

x  =  0.353  microns  (measured  in  units  of  '  at  1.G6  microns) 

(ThFu  film;  n  »  1.3-i8) : - ,  experimental  data:  - , 

theoretical  curve- 


Danage  threshold  data  and  theory  vs.  film  thickness  -  II] 

X  =  1.06  microns  (measured  in  units  of  \  at  1.06  microns) 

(ThF4  film;  n  =  1.3-il3) : - ,  experimental  data;  - 

theoretical  curve. 
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square  root  of  the  pulse  length,  but  that  sob*  other  material-dependent 
properties  enter  that  are  not  strictly  thermal  (because  of  the  wave¬ 
length  dependence  of  the  constant) .  In  other  words,  the  raacfaanisns  for 
absorption  in  films  of  various  fluorides  are  sufficiently  similar  that 
the  thermal  properties  dominate  the  damage  evolution,  while  oxides 
differ  from  film  to  film  in  some  way  other  than  the  obvious  thermal 
properties.  This  sane  behavior  is  also  exhibited  in  Fig.  21  for 
fluoride  and  oxide  data  from  a  different  source.  Note  hew  the 
fluorides  group  about  a  single  line  while  the  oxides  spread.  Mote  also 
that  the  order  in  which  the  oxide  lines  fall  is  the  same  in  both  data 
sets. 

When  the  fluoride  data  are  plotted  against  the  scaling'  (Eq. 

(10))  for  a  given  wavelength,  films  that,  are  full-wave  (at  1.06  ^m) 
have  damage  thresholds  that  are  similar  to  those  of  half-wave  films, 
but  as  we  consider  quarter-wave  or  snail  films  the  difference  becomes 
greater.  This  difference  may  be  explained  by  reference  to  Figs.  17-19. 
These  are  plots  of  the  exact  solution  of  the  diffusion  equation  using 
Hie  scattering  theory  to  compute  the  absorption  cross  section  of  these 
highly  absorbing  regions.  That  is,  a  combination  of  electromagnetic 
theory  and  thermal  diffusion  nicely  explains  the  increase  in  damage 
threshold  for  thinner  films  at  fixed  incident  wave  lengths.  This 
behavior  may  be  simply  understood.  As  the  film,  and  therefore  the 
ihsorbirc  regions,  becomes  small  with  respect  to  the  wavelength,  the 
regions  absorb  proportional  to  their  volunes.  However,  regions  of  the 
order  of  a  bciLf  wavelength  or  greater  absorb  in  proportion  to  their 
geccsetr  ilsii  cross  section  as  shown  in  Chapter  2.  Thus  as  the  film 
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Fig.  21.  Scaling  of  the  exper imenta  1  ly  measured  daoage  threshold  of 
both  fluorides  and  oxides  vs.  a  theoretically  derived 
parameter . 

The  parameter  contains  the  thermal  properties  and  tne  data  is 
for  half-wave  filss  at  *  =  0.248  sucrons  (*  measurements 
taken  at  20  ns) . 
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thickness  decreases,  absorption  decreases  with  the  geometrical  cross 
section  of  the  absorbing  region  up  to  seme  point.  Farther  decrease  in 
the  film  thickness  forces  the  absorption  to  decrease  more  like  the 
volune,  yielding  higher  damage  thresholds. 

'fhis  provides  the  relationship  between  the  size  of  the 
absorbing  region  and  the  wavelength.  In  addition  to  that  there  is  an 
increase  in  damage  threshold  due  to  the  fact  that  the  film  becomes 
thinner  than  the  diameter  of  the  absorbing  region  of  a  size  most  likely 
to  be  damaged.  Regions  smaller  than  this  require  a  larger  amount  of 
incident  energy  to  be  damaged. 

.  The  different  behavior  of  oxides  versus  fluorides  suggests  at 
least  one  of  three  things.  Either  there  is  a  major  oxide-dependent 
difference  in  thin  film  properties  as  opposed  to  bulk  properties,  or 
there  are  structural  or  chemical  complexities  not  accounted  for  or 
there  is  singly  a  mechanism  for  absorption  in  oxides  such  that  the 
absorption  appears  to  be  highly  material  dependent.  The  last  two  seem 
more  likely  because  of  the  wavelength  dependence  for  the  missing 
material  constant. ' 

It  should  be  pointed  cut  here  that  the  sore  predictable 
behavior  of  the  fluorides  does  not  necessarily  prove  that  an  absorbing 
inclusion  exists.  However  it  strongly  suggests  a  central  absorbing 
region  cocaon  to  the  various  fluorides  that  thermally  diffuses  the 
absorbed  energy  such  that  damage  thresholds  for  fluorides  are  thermally 
dominated.  This  region  should  have  properties  that  either  are 
initially  different  from  the  curroijnding  film  or  are  quickly  altered  by 
an  enhanced  fieio.  a  mechanical  or  chemical  process,  or  an  increased 


tenperaturg  in  the  region  or  possibly  sane  combination  of  these.  An 
enhanced  field  could  be  caused  by  a  microcrack  or.  a  local  defect 
[Bloenbergen  1973]. 


The  theoretical  evidence  that  the  size  of  the  absorbing  region 
that  is  danaged  is  of  the  order  of  the  film  thickness  does  not  exclude 
the  possibility  that  a  smaller  absorbing  region  expands  to  a  size  cf 
the  order  of  the  thickness  of  the  film  as  damage  ensues.  An  expanding 
absorbing  region  has  been  suggested  by  Kcoolov  [1982]  and  Anisimov  and 
Makshantsev  [1973]  and  expanding  damage  regions  have  been  observes  by 
Fcltyn  and  Jolin  [1933].  Also  Babadzhan,  et.  al  [1982]  have  suggested 
a  mechanism  for  this  process,  although  their  assumptions  leading  to  the 
absorption  cross  section  appear  to  be  in  error. 

It  is  interesting  to  note  the  stray  ThF^  points  in  Fig.  14  of 
the  quarter-wave  fils-.  These  same  points  are  plotted  on  Fig.  19  for  a 
damage  threshold  at  X  *  1.06  micros  versus  film  thickness.  They  are 
anomalously  high  on  this  plot  also,  indicating  that  they  are  unique 
points  which  do  not  follow  the  trends  of  the  surrounding  data. 

It  should  be  pointed  out  here  that  there  are  fundamental 
differences  between  oxides  and  fluorides  which  exhibit  themselves 
through  the  macroscopic  properties  of  the  materials.  For  example ,  the 
index  of  rifracticn  of  fluoride  films  is  usually  less  than  that  of 
oxide  films.  T he  melting  temperature  is  also  less  for  fluorides  than 
for  oxides.  Furthermore  the  UV  cut-off  of  absorption  occurs  for 
shorter  -wavelengths  iu  fluorides  than  in  oxides  and  thus  multi  photon 
effects  would  appear  sooner  in  oxides  as  the  wavelength  decreases. 

There  are  other  possible  differences  such  as  the  structural  properties. 
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the  elasticities  or  the  theraanecbanical  response  of  the  materials. 

One  final  apparent  difference  is  the  fact  that  fluorides  probably  have 
oxide  securities,  but  the  nature  of  the  iascurities  in  oxide  films  is 
unknown.  The  implication  or  significance  of  these  various  differences 
is  unknown.  However,  the  fact  that  oxides  tend  to  have  lower  damage 
thresholds  yet  higher  melting  tsseeratures  indicates  that  setae 
additional  mechanists  say  intervene  in  the  dsss^e  evolution  process  for 
oxides. 

One  additional  comparison  of  the  thenssd  scaling  is  given 
pjsrely  Her  sake  of  cccpletscess  (Fig.  22}.  This  is  plotted  along 
with  data  fros  {Rainer,  et  al.  1582}  for  toe  porpose  of  ccoxsrison. 

Sote  that  the  data  is  plotted  '  *■  different  scales.  Comparison  of 

these  scales  for  wavelength  variation  cay  or  oay  not  be  meaningful 
since  they  are  from  two  different  experiments. 

The  final  experimental  caparison  given  is  that  of  the  work 
presented  in  Chapter  3.  Sosa  preiisucary  results  fran  experiments 
currently  being  conducted  fey  Kardach  [1985}  shows  a  scaling  with  the 
themed,  conductivity  of  the  substrate  for  single  layer  films.  The 
results,  shown  in  Fig.  23,  are  consistent  with  the  findings  in 
Chapter  3. 

It  is  felt  that  significant  correlation  has  been  shewn  with 
experiment  of  the  concept:  dssaispad  ha re.  transport  scalings  in 
Chapter  6  indicate  that  the  absorption  oast  be  strong  ana  localized. 
This  is  further  si  reported  by  the  morphology  observed  in  tie  damage 
data.  The  next  chapter  concentrates  on  the  development  of  the  physical 
processes  leading  to  the  strong  localized  absorption. 


pig.  22.  Scaling  of  the  experimentally  measured  damage  threshold  of 
fluorides  vs.  a  theoretically  derived  parameter  -  IV. 

the  parameter  contains  the  thermal  properties  and  the  data 
[Bettis,  1975]  is  for  haif-*ave  films  at  1.06  microns  and 
40  ns.  Set**  of  the  uncoaaon  fluorides  in  the  data  set  did 
not  fit  the  scaling. 
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Fie.  23.  Dasage  thresbo3xi  vs.  substrate  thermal  conductivity. 

The  data  was  taken  at  1.06  microns  and  5  ns  on  half  wave 
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mSMSHTAL  ABSORPTION  HBCHAKISMS 

In  the  work  set  forth  to  this  point,  the  concepts  of  a 
localized  absorbing  region  within  an  optical  thin  film  have  been 
presented,  these  localized  absorbing  regions  are  assumed  to  be  due  to 
anomalies  of  so me  sort.  Heat  transfer  equations  were  solved  indicating 
trends  in  damage  thresholds  as  a  function  of  material  and  experimental 
parameters.  Thermal  diffusion  was  assumed  as  the  dominant  controlling 
mechanism  in  the  damage  process.  Experimental  data  correlates  quite 
well  with  this  model.  The  thermal  diffusion  equations  incorporate  a 
ncnhomogenecus  thermal  source,  which  is  due  to  a  localized  electron 
density  interacting  with  both  the  incident  field  and  the  lattice 
structure. 

In  actuality,  there  are  tao  coupled  nonfeomogeneous  diffusion 
processes  involved.  They  are  thermal  diffusion  due  to  t be  afore¬ 
mentioned  thermal  source,  and  electron  diffusion  due  to  a  localized 
electron  source.  The  electron  source  properties  depend,  in  the  train, 
upon  the  local  optical  field  and  local  temperature.  These  coupled 
equations  are  generally  (intractable  and  as  such  are  usually  decoupled 
through  assumptions  which  can  apply  in  certain  limits. 

Assumptions  that  were  made  are  that  the  operative  coefficients 
and  the  source  for  thermal  diffusion  are  independent  of  temperature, 
and  the  electron  density  has  been  assured  constant.  Thar  is,  ’-he 
electron  density  is  modeled  as  a  step  function  .,1  time  and  space  seen 


feat  it  evolves  instartaneGusly  and  does  Dot  diffuse.  These 
assumptions  are  obviously  not  valid  in  all  cases  of  interest.  Certain 
refractory  oxide  film  damage  [Walker ,  Guenther  and  Sielsen  1981} ,  for 
example,  exhibits  morphology,  which  grows  spatially  during  irradiation 
to  a  size  comparable  to  the  incident  beam  diameter. 

In  addition  to  the  short  casings  of  these  employed  assumptions, 
the  use  of  an  instantaneous  electron  density  must  be  justified.  Thus, 
mechanisms  for  fee  evolution  of  an  electron  density  have  to  be 
investigated . 

Theory 

-  c~ 

In  an  attempt  to  refine  the  model  of  Chapters  2  through  4  in 
the  direction  of  acre  cceprebensive  theory,  fee  concentration  of  effort 
is  placed  upon  the  evolution  of  the  localized  electron  density  and  the 
absorption  processes  involved.  The  approach  taken  to  this  end  is  to 
write  down  the  classical  equation  of  electron  diffusion  and  solve  it 
for  a  generalized  electron  source.  Conservation  of  electromagnetic  and 
thermal  energy  is  then  employed  wife  the  principal  purpose  of  ascer¬ 
taining  useful  scaling  laws.  The  comparison  of  these  scaling  laws  with 
experimental  data  implies  certain  preferred  absorption  mechanises. 

These  mechanisms  are  then  further  investigated  via  the  appropriate 
expressions  governing  electron  diffusion. 

The  principal  findings  are  feat;  i)  local  metallic 
concentrations  suppleuented  by  color  centers  may  be  a  viable  mechanism 
for  the  initial  electron  production  in  certain  optical  thin  films  and 
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ii)  very  sail  regions  of  these  anomalies  (  200  A)  can  possibly  lead  to 
danage  when  supplemented  by  an  avalanche  ionization.  Avalanche 
ionization  is  only  assured  to  occur  in  optical  coatings  and  serves  to 
expand  the  region  that  is  absorbing  at  the  anomalously  high  rate. 

The  classical  equation  of  electron  diffusion  is  given  asi 


'n  l  32 

j§-  f -g-j  UDen>  +  g{n,T,j£J)  -  (n  -  n.)Y(n,T> 


where  Oft  *  electron  diffusivity  (cm  /S),  |Ej  *  electric  field  strength 
(V/ca) ,  g  »  electron  carrier  generation  (cm  ^s” 1 ) ,  n  =  electron 
carrier,  density  (aa~^) ,  ni  ■  initial  carrier  density  (ci*5) , 

T  »  temperature  (  K) ,  and  Y(nfT)  *  free  carrier  decay  rate  (s  *) . 

The  assumptions  made  in  solution  of  the  above  equation  are  that 
De  is  independent  of  n  and  T  (which  is  not  true  in  general) ,  and  g  is 
independent  of  T.  the  dependence  of  g  on  n  and  E  are  incorporated 
iteratively  and  the  other  assumptions  will  be  addressed  in  Chapter  7. 

The  solution  of  the  above  equation  in  all  space  is  found  via 
integral  transform,  is  absorbed  in  n  so  that  n  becomes  n  *  n  -  n. . 
The  general  solution  is  found  to  be 


,  r-  -CDe52  >  Y)t 

n(r,t)  =  —  I  ds  sin(sr)e 
J 
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rz  r  tDe 

dt1  i  r'dr*  sin(sr‘)g(r' ,t  )e 
•<n  Jn 


(D i  S*+r)t' 


In  order  to  decide  upon  a  reasonable  choice  for  the  electron 
source  generation  tern  g  (r,  t) ,  attention  is  diverted  to  a  consideration 
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of  tha  conservation  of  energy.  The  morphology  of  damage  suggests  a 
localised  source  of  electron  generation,  however,  it  does  net  imply  a 
mechanics.  Scaling  laws  derived  from  the  conservation  of  energy 
consideration  nay,  though. 

The  Poynting  theorem  may  be  written  as 


|“=,7.s-(j  +  4!>-£4-S.i£  us) 

rtbere  w  is  the  field  energy,  S  is  the  Poynting  vector,  J  is  the 
electron  current  flux,  P  is  the  material  polarization  and  M  is  the 
magnetization  vector.  If  the  field  is  in  approximate  steady  state  with 
the  absorption  process  and  magnetic  transitions  are  neglected,  then 


^  ^  ^  3P 

»-5.J-B4S*{§.  (49 

The  J  •  £  term  contains  the  energy  stored  in  electron 

acceleration  along  with  electron-phonon  collisions  and  the  P  •  E  terra 

contains  the  energy  stored  in  polarization  along  with  electronic 

transitions.  The  J  •  £  term  is  the  dominant  term  in  avalanche- 

« 

ionization  and  inverse  brehrastrahlung  processes  while  the  P  -  1  terra 
is  assumed  dominant  in  single  and  aultiphoton  transitions.  All  of 
these  are  possible  mechanisms  of  absorption.  In  order  to  avoid 


restriction  to  any  certain  mechanism  the  v  •  §  term  is  used  to  represent 


absorption. 
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For  simplicity  an  incident  plane  wave  is  assumed,  which  gives 
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where  I,  n  and  u  are  the  incident  field  intensity,  real  index  of 
refraction  and  magnetic  permeability,  a  is  the  generalized  absorption 
coefficient  and  c  is  the  generalized  absorption  cross-section  with  V 
the  volume  that  is  absorbing. 

From  the  first  law  of  thermodynamics  the  energy  balance  of  an 
absorbing  region  is 

7  •  (-  7(k?))  +  (ocvT)  *  -  *  '  S  ,  (51) 

where  k  *  thermal  conductivity,  cy  *  specific  heat  at  constant  volisae, 
and  o  *  mass  density,  which  leads  to 

EjjCJ/ca2)  *^“J  %  (pcvT)  -  v  •  (v  (kT) ) }  ,  (52) 

where  is  the  energy  density  of  the  incident  pulse  at  which  the 

region  reaches  a  critical  temperature  which  causes  damage. 

The  energy  depletion  from  the  incident  field  occurs  over  a 

region  represented  by  the  length  a-1  which  includes  absorption  by 

electrons  that  ironed  lately  thensalize  along  with  those  that  thermal  ize 

after  electron  diffusion  of  length  6  .  As  oointed  out  by  Meyer, 

e 

Bartoli  and  Kruer  (1980} ,  the  thermal  diffusion  may  be  incorporated  in 
the  same  manner.  That  is,  the  energy  that  is  deposited  in  the  lattice 
is  deposited  in  a  region  resulting  from  absorption  and  thermal ization 
of  electrons  (both  with  and  without  diffusion)  along  with  the  thermal 
diffusion  6^.  that  occurs  after  deposition  during  the  pulse.  The 
thermal  diffusion  length  is  given  by  =  »t)t  .  The  rwo  parallel 
processes  may  be  added  to  give 
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(53) 


Here  Xj,  is  the  fraction  of  laser  energy  deposited  directly  into 
the  Lattice  and  x9  As  the  fraction  of  energy  that  diffuses  before  being 
deposited. 

•Thus  we  can  write 
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by  substituting  an  effective  diffusion  length  for  the  diffusion  terra. 

The  importance  of  this  procedure  is  the  dernonstratioju-of  useful 

scaling  laws.  For  exs^le,  when  6_  >  d  ,  a_i, 
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The  first  case  exhibits  idle  same  scaling  law  derived  for  a 
spherical  absorbing  region  where  thermal  diffusion  dominates  (see 
Chapter  2).  Here,  it  is  much  more  general  in  that  no  specific  geometry 
has  been  assumed  and  temperature  variation  of  the  constants  is 
incorporated. 


The  second  case  demonstrates  an  absorption  dominate  1  process 
and  implies  that  <=  a”1.  Thus  the  wavelength  dependence  o>.'  a  (A) 
should  exhibit  itself  in  the  damage  threshold  data  roughly  as 
Eq  «  a"1  (A). 

In  addition  to  the  exhibition  of  wave  length  dependence  in 
certain  cases,  an  important  point  of  the  above  scalings  is  the 
demonstration  of  the  perspective  which  tl  •'  lend  to  the  all  important 
energy  transport  process.  The  scales  over  which  energy  transformation 
and  transport  occur  are  solely  responsible  for  whether  the  damage  is 
absorption  or  diffusion  dominated. 

The  ’wavelength  dependence  for  damage  from  reference  [Walker,  et 

al.  1981]  may  be  found  by  averaging  the  wavelength  dependence  of  the 

films  (A  and  A/2)  for  five  and  fifteen  nanoseconds.  Wien  this  is  dons 

it  is  found  that  between  A  *  0.26  microns  and  0.53  microns  the 

1  95  2  49 

fluorides  damage  as  A  *  while  oxides  damage  as  A  ,  ana  between 

0  24 

A  *  0.53  microns  and  1.06  ah  crons  fluorides  damage  as  a  *  and  oxides 
A0’5.  There  is  significant  statistical  deviation  frcm  film  to  film 
(i.e,  for  the  various  coating  materials). 

In  all  of  the  theories  investigated,  tnis  increasing  power  of 
frequency  in  the  absorption  furcton  with  increasing  photon  energy  was 
only  noted  in  photo- ionization  processes.  However,  the  films  and 
photon  energies  investigated  [Walker ,  et  al. ,1981]  should  not  be  near 
band  edge  transitions.  It  should  be  noted  that  a  large  density  of 
impurity  states  can  exhibit  this  sane  behavior.  In  addition,  the 
transition  between  different,  but  related,  processes  of  absorption  at 
different  frequencies  is  a  possibility  (i.e.,  different  sorts  of 


impurities) .  These  may  occur  in  the  same,  or  even  separate  discrete 


regions. 

The  oeneral  a  is  constituted  as  a  =  )  a  +  ac  where  7  a,  is 

*■  n  Fc  -  n 

n  n 

the  absorption  coefficient  for  n-pboton  absorption  and  apc  is  the 

absorption  coefficient  for  free  carriers.  £  a  is  related  to  apc  ir 

n 

that  n-phcton  absorption  produces  free  carriers. 


Although  these  processes  are  generally  considered  intrinsic 
they  also  may  apply  to  the  obviously  extrinsic  damage  processes  which 
occur  in  optical  coatings  due  to  localized  shallow  or  impurity  states 
in  the  coating.  These  shallow  impurity  states  support  photo-ionization 
and  impact  ionization  at  threshc—  levels  far  belcw  those  expected  of  a 
pure  material.  These  states  may  occur  due  to  imparity  atoms, -dons, 
non-.sto ichicxnetr ic  mixtures  leading  to  metallic  colloids,  color  centers 
qc  simply  because  of  local  disorder.  Shallow  states  due  to  disorder 
5 i«ve  been  confirmed  for  D.C.  fields  by  conductivity  measurements 
[Ovshinsky  1563 3 

Thus,  due  to  the  above  possibilities  and  the  observed 
morphology  it  is  assuaed  that  there  exists  the  probability  cf  very 
small  localized  regions  of  impurities. 

One  obvious  choice  for  an  electron  source  function  is  a  uniform 

spherical  region,  of  impurities  or  shallow  states,  of  radius  r  .  The 

0 

origin  of  su-J.  a  region  say  be  a  coagulation  of  centers  resulting  from 
dislocations,  as  discussed  below. 

This  source  term  is  given  by 


Where  for  single  and  multi photon  transitions  of  shallow  states 


fit'}  =* 


1  hv 


.Ik, 

2hv; 


•3 (t* ) • 


(58) 


Here  I  is  the  incident  intensity  (ev/c a  -s) ,  is  the  one  photon 

_ eoeff ieient  (l/cm) ,  s  is  the  two  photon  coefficient  (cs-s/ev) ,  ip  is 

the  quanta  efficiency  (1),  hv  is  the  photon  energy  (ev) ,  and  the  ^’s 
are  steo  functions.  This  integral  (Eg.  (57)  into  Sg.  (47/)  is  very 
similar  to  the  one  in  Chapter  2  and  affords  a  solution  of 
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cc-  the  region  of  parameter  space  of  interest  here,  the  above 
integral  provides  raxaerical  problems  in  its  evaluation.  However  the 
first  j-p?™  of  the  integral  nay  be  done  analytically  giving  a  solution 

of 
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In  alkali -tel  ides  the  shall  a*  inpurity  states  may  be  provided 
by  a  distribution  of  metallic  colloids  res’ll  ting  from  poor  stoichio¬ 
metry  and  various  aggregates  of  color  centers.  Colloids  may  be 
produced  in  tie  process  of  heating  to  200  to  350  C  [Schulman  and 
Compton  1963;  and  Hayes  19741  even  by  insufficient  heating  during 
film  deposition,  in  addition,  mild  beating  causes  the  coagulation  of 
F-centars  to  colloidal  particles  [Schulman  and  Ccscton  1963] .  With 
specific  reference  to  alkaline  earth  fluorides,  HF-centres  in 
addi tively-color ed  alkaline  earth  fluoride  crystals  readily  aggregate 
forming  more  carplex  structures"  [Hayes  19741.  An  example  of  the  sort 
of  distribution  that  say  occur  is  shewn  in  Fig.  24, 

Where,  in  optical  coatings,  the  possible  localized 
concentrations  of  ai  xirpticn  centers  such  as  these  cone  from  is  not 
discussed  here.  Dislocation  clusters  and  poor  stoichiometry  are  quite 
likely  in  optical  coatings.  The  assumption  made  is  that  the  conditions 
of  deposition  and  processing  of  coatings  allow  this  to  be  true. 

As  an  illustrative  example  of  one  of  the  materials  of  interest 
[Walker,  et  ai.  1981] ,  Ca?2  has  an  absorption  peak  at  521  no  in  what  is 
called  the  S  band  predominantly  due  to  fi-centers.  It  should  be  pointed 
out  that  in  CaF2  the  F-centsrs  peak  at  376  na  and  certain  colloids  peak 
at  550  m.  In  addition  to  these  are  many  other  complex  combinations  of 
centers.  [Scfeuinan  and  Compton  1963;  Hayes  1974;  and  Fowler  1568] . 

The  oscillator  strength  nay  be  assumed  to  be  %  1  [3eatracnt  and 
Hayes  1969;  and  Fowler  1564'  and  the  transition  rate  for  single  photon 
transitions  is  given  as  [Yariv  l'.'~  ^ay  L?73aj 


COLLOID  ABSORPTION  BANDS 
U.»FtS+.  R  ECT  CENTERS 
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Fig.  24.  A  possible  absorption  distribution  vs.  photon  energy 
for  a  couples:  aggregate  of  absorption  centers. 
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feik  3  ochv  ' 

where  g(v)  is  the  normalized  lineshape,  f ■ 


(64) 

is  the  oscillator  strength 


sad  1  is  tie  laser  intensity.  g(v)  [Ysriv  1967;  and  Gray  1973b]  can  be 
expressed  as  follow. 


g(y) 
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One  of  tise  wavelengths  of  interest  is  530  n o  [Walker-  et  al. 
1331] ,  which  is  slightly  off  center  of  the  transition  line.  Due  to  the 
large  width  of  the  lines  [Schulcaaa  and  Caspton  1963;  Kayes  lc,4;  Fowler 
1968;  Gray  1973b;  Sparks  and  Duthler  1975],  we  find  g(v)  =  0.9 ft", 
however  to  be  conservative  a  line  shape  factor  of  0.1/iv  will  be  used. 

Elan  this  arputatico  It  is  JoraJ  that  «ik  -  l.WV1  £tan 
ground  state  to  the  first  excited  state.  This  first  excited  state  is 
long  lived  [Beauaont  and  Hayes  1969;  and  Fowler  1964]  and  will  easily 
ionize  deeply  into  the  conduction  band  frcsa  here. 


Kesults 

Assicing  in  absorption  cater  density  of  a  3  10* 7  na'-5  (which 

21 

nay  be  conservative  for  a  thin  fiia  [Sparks  and  Duthler  1975]) 


w 

hv 


—=•  *  7x10 


(65) 


Haze  a  factor  of  1/2  was  incli  to  account  for  the  final  i  .ization 
step,  souse  o£  which  say  occur  thermally.  This  factor  say  as  placed 
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into  the  quantun  efficiency  n_  •*  1/2  and  can  be  considered  conservative 
due  to  the  ionization  transition  rate  approximated  from  equations  in 
references  {iariv  1967 ;  Vaidyanathan,  et  ai.  1977;  Vaidyanathan ,  et  al. 
1579]. 


i 

i 

f 

i 

i 

i 

i 

[ 


Several  important  points  related  to  Table  4  be 

mentioned.  The  results  that  follow  are  based  upon  classical  diffusion 
which  should  n ere  apply  where  stall  nanbers  of  electrons  are  involved. 

If  one  rail  tipi  ies  tie  electron  densities  by  the  imparity  region  volume 
it  is  clear  that  rather  small  numbers  of  electrons  are  indeed  involved. 
There  is  no  attempt  made  here  to  model  this  situation  rigorously.  This 
aspect  of  the  study  is  used  to  demonstrate  relative  orders  of 
magnitudes,  and  the  feasibility  of  very  small  local  anomaly  regions 
producing  large  electron  densities  in  the  conduction  band. 


Some  of  tbs  numerical  results  from  the  integration  of  Sq-  (59) 
are  she*®  in  Table  4. 

a  *  4.91  (csT1) ,  (for  n  *  10A7csf3,  c  *  2.1xlQ~i7) 

3  “0. 

\  *  1/2 

*  0.01  (cs//s) 

f,,  -*  l.Q  (oscillator  strength) 
gfv)  «  0.1/Av  (line  shape  factor) 
hv  a  2.34  eV 

,  -3 

to  =10  s 


radius  of  region  of  impurities 


*  conduction  electron  recombination  rate  (s~^) 

2 

*  absorption  cress  section  (cm  ) 

»  probability  of  electron  production  within 

volume  of  interest 

3  electron  density  at  0.1ns  and  r  =  50  A 

*  time  at  which  probability  of  electron 

production  exceeds  1 
For  sxanple,  the  initial  color  center  density  for  the  table 

below  is  chosen  to  be  10i/  era”3,  however  the  final  electron  densities 
produced  free:  this  exceed  this  number.  This  number  is  not  chosen  to  be 
realistic  as  ouch  as  it  is  chosen  to  be  conservative.  In  a  region 
where  trap  centers  ate  generated  by  polishing  damage  or  sane  other 
process,  the  trap  canters  would  probably  be  orders  of  magnitude  greater 
than  10*  esa  ,  Thus,  IQ  is  chosen  to  provide  a  conservative  transi¬ 
tion  rate,  though  it  is  not  allowed  to  limit  tbs  nuscer  of  electrons 
which  transition.  Conduction  electrons  are  limited,-  however,  by 
recasbinaticn,  or  decay.  It  may  further  be  noted  that  if  IQ17  cnf J 
color  centers  existed  in  a  region  of  100  A  in  radius,  at  best  one  half 
of  an  available  electron  would  be  in  this  region.  Clearly  its  only 
purpose  is  to  give  conservative  transition  rates.  As  such,  the  numbers 
in  'Sables  4  and  5  should  not  be  taKen  as  accurate.  They  are  presented 
to  demonstrate  relative  orders  of  magnitudes  of  conduction  electrons 
that  could  be  produced  under  the  conditions  of  radiation,  decay  ana 
diffusion.  One  prisary  point  of  Table  4  is  the  demonstration  of  the 
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Table  5.  Electron  Densities  achieved  at  12«rI0'^  cm  from  color  center 
initiation  dm  to  aa  electron  avalanche  for  various 
parameters.  The  avalanche  is  assumed  to  spread  via 
diffusion. 


=  lO10^"1} 


De  =*  0.01  (cm  /s) 


£i* 


1.0 

g(v)  *  0.1/Av 
hv  *  2.34 


10 


10 


20 


Avalanche  coefficient 
Electron  diffusion  coefficient 
Oscillator  strength 
Line  shape  factor 
Photon  energy 
Pulse  length 

Limiting  electron  density 


*  The  limiting  electron  density  n  was  chosen  as  a  reasonable  value 
based  upon  [Burs te in  1954 J. 
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fact  that  the  size  of  a  region  and  the  intensity  of  radiation  dictate  a 
threshold  condition  due  to  carpeting  probabilities  of  transition  and 
decay  in  the  region.  'What  constitutes  this  region  is  cot  assumed  here, 
other  than  the  fact  that  it  contains  shallow  centers.  It  may  be  an 
inclusion  at  the  base  of  a  nodule  [Lao  and  Saith  1935]  or  at  a  grain 
boundary.  It  nay  be  a  dislocation  cluster  [Lewis,  et  al.  1985].  It 
should  be  noted  here  that  colloidal  particles  deposit  selectively  at 
dislocations  [Schulman  and  Compton  1363} .  In  addition,  dislocations 
are  the  source  of  color  centers. 

This  threshold  is  not,  however  the  damage  threshold.  The 
important  parameter  in  damage  is  the  rate  at  which  lattice  energy 
builds  in  a  local  region.  Even  if  the  electron  transition  threshold  is 
exceeded,  the  energy  say  be  harmlessly  conducted  away  by  electron  and 
phonon  transport.  Table  4  also  serves  to  denonstrate  that  if  suffi¬ 
cient  nunbers  of  electron  trap  centers  (M-centers,  in  this  case)  are 
present,  the  transition  is  quick  («  t  ) .  The  partitioning  of  such 

p 

regions  where  the  volume  say  be  smaller  than  the  volume  between  trap 
centers  is  very  artificial.  It  is  intended  to  generalize  the  table  to 
weaker  transitions  at  higher  densities,  or  perhaps  colloidal  transi¬ 
tions.  This  is  why  the  cross  section  is  listed  in  Table  4.  Whatever 
the  source  of  electrons,  if  the  cross  section  is  as  given,  the  table 
should  remain  valid  in  the  region  of  the  transitions. 

Although  the  electron  densities  represented  in  these  small 
regions  are  quite  high,  they  could  not  alone  transfer  sufficient  energy 
to  the  lattice  to  cause  damage.  Densities  of  this  sort  are  required 
over  regions  approaching  0<1  to  i.  micron  in  diameter  [see  Chapter 
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Two] .  Since  anomaly  regions  of  this  size  are  not  observed  in  optical 
coatings  (except  for  nodules} ,  these  absorbing  regions  must  evolve 
during  irradiation.  Since  nodules  are  of  the  sacs  material  as  the 
surrounding  film  the  absorption  in  and  around  them  must  similarly 
evolve.  Of  the  various  mechanises  studied  for  this  process,  the  only 
reasonable  explanation  found  is  the  expansion  of  the  absorbing  region 
via  impact  ionization.  That  is,  these  snail  local  high  densities  of 
conduction  electrons  ionize  the.  surrounding  off  resonant  electron  trap 
centers  and  impurities.  This  is  casaoniy  termed  avalanche  ionization, 
though  it  usually  applies  to  purely  intrinsic  processes  (ionizing 
valence  electrons} .  It  should  be  noted  that  if  the  absorption  center 
region  is  on  the  order  of  0.1  to  1.  micron,  no  avalanche  is  required  to 
produce  the  appropriate  region  of  high  electron  densities.  This, 
however,  is  probably  unlikely.  A  solution  of  this  case  is  shown  in 
Fig.  25  which  emphasizes  the  effect  of  y. 

An  avalanche  due  to  impact  ionization  of  impurity,  colloidal  or 
trap  states  would  explain  the  observation  of  lower  dsaage  thresholds  in 
doped  materials  regardless  of  whether  or  not  the  doping  impurities  are 
in  resonance  with  the  field.  The  pbotoionization  of  centers  to 
initiate  the  process  «3uld  explain  qualitatively  the  wavelength 
dependence  of  damage  thresholds  in  optical  coatings. 

Although  specific  nunbers  are  used  for  a  given  case  in  a 
certain  material,  this  is  done  for  cccputational  purposes.  It  is  felt 
that  the  calculation  of  a  specific  exacple  is  desirable.  The  concepts 
are  intended  to  be  applied  to  more  general  cases  through  the  use  of 
absorption  cross  sections. 


A  detailed  study  of  avalanche  ionization  will  not  be  conducted 
hare.  That  has  been  the  topic  of  considerable  research  by  authors  such 
as  [Sparks  and  Duthler  1975;  Yablonovitch  and  Bloeraberqen  1972;  Fradin, 
Yablcoovitcfa  and  Bass  1973;  Epifacov  1974] .  As  reference  [Fradin,  et 
al.  1973]  implies,  avalanche- ionization  has  been  confirmed  in  solids, 
and  as  reference  [Sparks  and  Duthler  1975]  indicates,  a  computational 
competency  has  been  achieved  with  the  processes  involved.  It  will 
simply  be  assumed  here  that  due  to  shallow  trap  states,  densities  in 
t-ccess  of  those  required  [Sparks  and  Duthler  1975}  to  initiate  an 
avalanche  have  been  provided.  An  exponential  growth  then  follows, 
which  is  subsequently  limited  by  ionizaticg  of  a  large  fraction  of  the 
available  atoms,  carrier  decay  and  the  dynamic  Burstein  shift  [Burstein 
1954] .  A  simple  calculation  with  Bq.  (67)  shows  that  at  optical  fre- 
quencies  (^500  rxn)  and  damage  intensities  (^1-2  Gw/cnO  in  fluoride 
films  a  conduction  electros  can  gain  energies  of  one  naif  to  three  sV 
in  1G~10  seconds. 


de 

dt 


* 

a 


(1+cA^2) 


(67) 


(In  Eq.  (67)  z  is  tbs  energy  transferred  to  the  elec-trcn  and  is  the 
time  between  collisions) . 

The  spread  in  values  is  due  to  uncertainties  in  collision 
frequencies*.  These  energies  are  sufficient  to  ionize  many  electron 
trap  centers  and  seme  impurities  [Sparks  and  Duthler  1977}. 

The  simple  model  used  herein  will  produce  a  conservative 


estimate  of  the  growth  of  an  electron  density  initiated  by  a  small 
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region  of  absorption  centers.  This  relates  to  the  findings  of  Chapter 
2  which  demonstrate  a  requirement  of  isolated  regions  on  the  order  of 

the  film  thickness  with  an  average  electron  density  on  the  order  of 

,rt18  -3 
10  cm  . 

For  this  case,  the  form  of  g(r,t)  chosen  is 


g(r\f)  *  f(t'}e~p‘*(t,)r'2 


(68) 


where  p(t’)  represents  a  spatial  expansion  of  the  electron  source  due 
to  the  avalanche  process.  It  is  assuaed  from  previous  studies  [Fradin, 
et  al.  1973]  that  an  avalanche  region’s  growth  is  nonlinear,  however  a 

conservative  estimate  of  the  region's  size  at  a  time  t  is  5  *  •'D  t. 

e  --  e 

That  is,  the  distance  electrons  have  travelled  by  diffusion  from  their 
initiation  point.  Thun;  pit)  *  l/*tTt  can  be  used  as  a  tried 
parameter. 

Subsequent  integration  of  Bq.  (47)  affords 


n(r-t) 
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0  CD  ( (t  -  t’>  +  — ~  )) 
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«p(-r(t-f))  exp^-fcT)  — '  L 


4Djp‘ 


(69) 


with 


f(t')  = 


a  t’ 

a 


(7C) 


where  a  is  the  avalanche  coefficient.  Equation  (70)  simply  represents 
a  taeans  of  generating  electrons  initially  at  an  exponential  rate.  This 
growth  is  terminated  at  sane  point  due  to  processes  such  as  the  dynamic 
aurstein  shift  and  ionization  of  a  large  fraction  of  the  available 
electrons, 

A  rigorous  justification  of  Bqs.  (69)  and  (70)  above  is  not 
possible.  These  expressions  are  simply  a  reasonable  choice  for  the 
purpose  of  a  Gedanken  study  of  a  postulated  process.  It  should  also  be 
noted  that  as  shown  in  Eraden,  at  al.  [1973]  avalanche  is  not  a 
symmetric  process.  It  in  fact  tends  to  propagate  in  the  direction  of 
the  beam  source  due  to  an  attenuation  of  the  field.  In  optical  thin 
films  this  propagation  is  modified  by  encounter  of  the  film  surface. 

Sane  of  the  numerical  results  are  presented  in  Figs.  26  and  27 
and  Table  4. 

Discussion 

The  numerical  results  presented  in  the  previous  section 
indicate  the  feasibility  of  microdamage  sites  initiated  by  a  localized 
absorbing  anaaally.  An  avalanche  of  electrons  can  prcduce  large, 
densities  of  electrons  expanding  cut  to  diameters  on  the  order  of  the 
film  thickness. 

It  should  be  noted  that  the  field  thresholds  are  lower  than 
those  reported  previously  for  avalanche  processes  in  similar  materials. 
However,  those  were  reported  in  bulk  solids.  The  starting  mechanism  m 
thin  films  should  provide  lower  thresholds  for  starting  the  process 
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The  growth  of  the  electron  density  -  II. 

The  growth  of  the  electron  density  at  a  radius  of  100&  for 
various  parameters.  The  parameters  chosen  are  considered 
reasonable.  Electron  densities  becaae  very  high  at  this 
radius  bosh  with  and  without  avalancne  ionization.  Dg  is 
assumed  to  be  0.01  ctr/s. 
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(colloids  and  aggregates  of  high  density  centers) .  The  large  density 

of  shallow  states  (centers)  which  nay  coagulate  about  colloids,  or 

coagulate  about  themselves,  should  provide  lower  levels  required  for  a 

sustained  avalanche  process  (i.e.,  the  ionization  of  the  off  resonant 

centers  requires  significantly  less  energy  than  valence  electrons)  - 

9  10 

The  avalanche  coefficients  (10  -10  )  were  chosen  oy  scaling 

data  frcm  Yablonovitch  and  Bloenbergen  [1972]  using  relations  of 
Epifanov  (1974}  for  2.0  to  2.5  ev  inpurity  states-  These  numbers  were 
further  supported  by  Bq.  (67) .  It  should  be  noted  here  that  5q.  (67) 
is  overly  simplistic  and  nay  not  be  sufficiently  conservative.  In 
addition  Epifanov  [1985]  believes  that  the  Coefficients  determined  in 
Yablcrwvitch  and  Bloenbergen  (1972]  nay  be  too  high. 

The  approach  taken  here  nay  not  be  considered  rigorous.  This 
is  intended  to  be  a  feasibility  study.  The  feasibility  of  this  process 
or  a  similar  one  is  considered  sufficiently  demonstrated.  The  exact 
processes  involved  are  still  considered  unknown. 

Conclusions 

It  has  been  shown  that  due  to  the  non-ideal  structure  of  an 
optical  thin  film,  large  local  electron  densities  say  be  generated  by 
non-stoichiaaetric  colloids  or  aggregates  of  color  centers.  These  are 
likely  to  exist  in  optical  systems  due  to  the  conditions  of  deposition, 
substrate  polishing  and  subsequent  treatment,  if  avalanche  ionization 
does  occur  in  these  films,  regions  as  small  as  100  A  in  radius  of 
colloids  and  centers  nay  initiate  the  process.  Otherwise  regions 
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3  4 

closer  to  the  thickness  of  the  film  10  >10  A  of  colloids  and  centers 
may  be  required. 

It  should  be  pointed  out  that  colloids  in  the  fluoride  films 
would  probably  be  oxidized  due  to  ewer  present  atmospheric  oxygen,  the 
process  is  thus  further  complicated.  Band  edge  transitions  may 
actually  be  occur ing  in  oxidized  metallic  colloids  in  fluoride  films. 

In  addition  to  this,  correct  absorption  cross  sections  of  metal 
colloids  in  the  optical  region  have  yet  to  be  computed.  Errors  or 
inappropriate  assumptions  have  been  found  in  those  reported  in  the 
literature. 

Measurements  of  colloidal  absorption  bands  indicate  that  it  is 
likely  that  colloidal  absorption  at  X  *  1.06  micron  is  a  strong 
competitor  for  the  two-photon  process  which  would'  be  required  for 
M-centers  in  CaF? .  The  two-photon  computation  produces  transition 
rates  significantly  less  than  the  one  photon  at  0.53  micron.  However, 
correct  colloidal  cross  sections  neod  to  be  computed  for  comparison. 
Also,  more  rigorous  avalandje-computatiocai  need  to  be  performed  for 
this  specialized  case,  if  possible.  It  nay  not  be  possible  due  to  the 
lack  of  knowledge  of  the  thin  film  structure,  in  addition, 
complications  such  as  oxidation  and  electron  density  variation  of 
and  y  need  to  be  accounted  for. 

Despite  the  details  that  have  been  neglected  it  is  believed 
that  feasibility  has  been  shown. 


CHAPTER  7 


TQ-PERATURE  DEPENDENT  EFFECTS 

Theoretically,  temperature  dependence  exhibits  itself  through 
several  mechanisms.  The  first  and  most  obvious  is  through  the  electron 
-distribution.  Increased  temperature  means  increased  phonon  activity 
and  increase  in  electrons  statistically  promoted  to  the  conduction 
band,  from  either  a  shallow  state  os  a  valence  band. 

The  second  mechanism  for  temperature  dependence  is  the  band 
gap.  The  band  gap  invariably  decreases  with  increased  temperature. 

A  third  affect  of  temperature  is  its  a&vcu.'t  r?  tm  material 
parameters .  That  is,  the  thermal  diffusion,  electron  diffusion, 
electron  recombination,  etc, 

Cos  final  role  that  temperature  has  is  the  secondary  role  that 
it  plays  xn  the  primary  processes,  such  as  avalanche  ionisation.  Here, 
the  bespjc.ittfxe,  as  alssve,  exhibits  itself  via  electron  phonon 
collisions.  The  high  energy  griooons  assist  in  promoting  electrons  up 
La  energy  in  the  ccoesKtiee  sand.  This  role  is  considered  of  minor 
is^crtaoce  since  it  ragnizes  a  quantitative  adjustment  of  the 
dominant  processes,  it  zs  very  important,  however  for  quantitative 
calculation  for  comparison  with  data. 

A  model  of  the  electron  distribution  in  ideal  solids  is  well 
known,  and  can  be  shown  to  be  [Beam  1965} 
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n  -  Sisal-  m'VV3/i  e'^.  (71 

6  h3/2  “  P 

In  this  relation,  ne  is  the  conduction  electron  density,  k  is  the 
boltsrenn  constant,  taj^  and  a*  are  the  effective  electron  and  hole 
masses  and  h  is  Planck's  constant.  This  is,  of  course,  directly 
related  to  the  free  carrier  absorption  coefficient  by  the  Drude  model 
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where  u  is  the  magnetic  permeability,  tk  is  the  momentum  exchange 
(collision)  time,  u  is  the  incident,  field  frequency  and  nr  and  n&  are 
the  real  index  and  electron  density. 

Kanolov  [1982]  finds  experimentally  that 
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a(T)  »  Oj  +  aQ  e 


where  for  Ti02  at  optical  frequencies 
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a(T)  =«  (8Q  -  180)  +  (4.x107-108)  e“V21rco°-kl7f600)/T(aj,-1) .  (74) 
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Since  the  band  gap  is  at  about  3eV  *  2kT  this  yields  a  value  of  T  =  u 
*  l.SxlO4.  With  a  melting  temperature  of  2133  K  this  yields  a  maximum 
value  at  melting  or  a(T)  *  180  +  10 8  exp  (-17,500/2133)  *  2.6xl04. 
Although  this  is  nearly  on  the  order  of  the  required  absorption  ('-10”) 
it  requires  reaching  the  melting  terrpsrature  to  achieve  this  value. 
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The’ value  of  aQ  =»  10a  is  also  about  two  orders  of  magnitude  greater 
than  theory  predicts.  However,  the  temperature  dependence  of  the  band 
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gap  will  contribute  to  this  process. 

The  band  gap  variation  may  be  included  as  Eg  *  EgQ  -  -fib  where 

a  relatively  large  value  of  y  frcm  literature  is  8  [Popov  and  Fedorov 

1983}.  To  be  liberal,  the  value  y  =  10  is  chosen.  Thus  E  =  3eV  - 

10(0.86x10  4)  (2133)  *  1.17  at  melting,  so  that  a(T)  =  1.7xl05cm  1 . 

Under  these  liberal  conditions  it  is  clear  that  temperature  dependence 

can  be  important.  However,  for  temperature  dependent  absorption  to 

-4 

dominate  requires  an  initial  absorbing  impurity  on  the  order  of  10  cm 
to  initialize  the  process  [Kanolov  1982].  Absorption  centers  in  cannon 
oxide  coating  materials  (M^OjSiO^Al^O,)  are  well  known  by  those  who  ' 
study  color  centers  [Schulman  and  Compton  1963] .  It  is  also  likely 
that  many  impurity  transitions  are  within  range  of  the  conduction  band 
in  oxides.  Thus  if  absorption  centers  of  sane  sort  start  an  isolated 
absorption  in  a  film  such  as  TiO^  temperature  dependent  absorption  may 
contribute  via  a  spread  by  an  absorption  wave.  This  is  an  alternative 
to  an  avalanche  ionization  process.  It  is  likely  that  if  conditions 
are  appropriate  for  an  avalanche,  it  will  dominate.  It  is  a  generally 
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faster  process  than  temperature  dependent  absorption  and  does  not  in 
theory  demand  such  a  significant  initial  absorption  required  to 
initiate  a  thermal  instability  involved  in  temperature  dependent 
absorption.  Furthermore,  temperature  dependent  absorption  ray  apply  to 
only  a  select  range  of  materials  with  sraaii  band  gap  [Kanoiov,  1982] . 
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For  example,  if  a  typical  fluoride  film  CaF?)  is  considered, 
it  can  be  found  that  a  =  102  +  108  exp(-5.8xi04/1775)  *  102+  8.35xl0~7. 

If  the  temperature  dependence  of  the  band  gap  is  included  E  *  lOeV  - 

-4  2  -=  -i  5 

10(0.86x10  )  .  (1775)  *  8.5,  a  »  10  +  8.1x10  Jan  ,  This  is  still  an 

insignificant  absorption  at  the  melting  temperature.  In  addition,  it 

would  require  a  very  significant  absorption  to  begin  with  in  order  to 

achieve  a  temperature  near  melting. 

Thus,  the  only  possible  primary  dependence  on  temperature  for 
this  case  would  be  in  the  assist  of  the  ionization  of  shallow  trap 
centers  and  impurities.  Since  photo ionization  of  resonant  trap  centers 
occurs  on  a  scale  of  10  11  seconds  it  is  not  likely  that  ionization  of 
these  is  thermal.  If  temperature  dependence  does  play  a  part 'in  the 
absorption  in  large  band  gap  coatings  it  is  more  likely  that  its 
contribution  is  in  the  assistance  of  impact  ionization  and  photon 
absorption  [Sparks  and  Euthler  1975 J.  Another  way  to  lock  at  this  is 
to  realize  that  for  high  intensity  optical  frequency  interaction  the 
electron  temperature  leads  the  phonon  temperature.  That  is,  the  field 
interacts  solely  with  electrons  which  impart  cheir  energy  to  phonons. 

It  is  less  likely  that  electrons  impart  their  energy  to  phonons  and 
equilibrate  with  other  phonons  and  then  ionize  valence  or  impurity 
bound  electrons  than  it  is  that  electrons  are  directly  ionized.  Trap 
centers,  by  their  nature,  have  larger  electron  cross  sections  than  the 
lattice  ions.  Any  significant  density  of  both  trap  centers  and  excited 
carriers  necessitates  impact  ionization.  Furthermore,  simple  computa¬ 
tions  (Eq.  (67) )  show  that  an  optical  field  (^  500  nm)  at  damage 
intensities  (^  1-2  (GW/an~))  for  fluorides  can  generate  one  half  to 
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three  eV  electrons  in  10  seconds  (depending  upon  the  coliisionai 
frequency) .  These  energies  are  sufficient  to  ionise  many  trap  centers 
1  to  3  eV)  along  with  sense  inpurities  (Sparks  and  Duthler  1977] .  In 
this  same  tine  period  even  che  most  strongly  absorbing  (a  *  105  cm  *) 
regions  could  increase  at  best  a  couple  of  tens  of  degrees  K  in 
temperature.  If  an  absorption  sue*  as  this  exists  there  is  no  need  for 
temperature  dependent  increase. 

The  temperature  dependence  of  the  material  parameters  nay  also 

1  r 

be  expressed  theoretically.  The  thermal  conductivity  k  *  —  I  c^L  , 

t 

where  i  denotes  the  type  of  carrier,  c  is  the  contribution  of  carriers 
to  the  specific  heat,  v  is  the  carrier  velocity  and  1  is  its  mean  free 
path.  At  low  temperatures  k  «  T~n  where  n  *  0,  1,  2',  or  3, "due  to  the 
phonoi  frequency  dependence  of  1  and  c  [Touloukain  1970] .  It  generally 
bas  a  positive  slope  at  lower  temperatures.  At  high  temperatures  k(T) 
is  primarily  governed  by  the  temperature  dependence  of  1  (i.e.  1  *1/T) . 

With  the  exceptions  of  S107 ,  HfO,  and  ZrO,  the  optical 
materials  (Walker,  1981]  are  on  the  down  slope  of  the  conductivity 
curve  at  roan  temperature  and  above,  with  usually  similar  slopes 
(Touloukain  1970] .  Hf02  and  Zr02  are  approximately  independent  of 
temperature  and  5i02  is  clearly  positive  {Z  Ti/ 2 ) .  Since  the  rata  at 
which  the  heat  is  conducted  away  by  the  surrounding  film  is  the  point 
of  interest,  and  in  the  problems  of  interest  significant  conduction 
takes  place,  it  is  felt  that  a  relative  measure  of  this  rate  is  given 
by  the  room  temperature  conductivity.  This  is  particularly  true  for 
materials  where  the  behavior  of  the  conductivity  with  temperature  is 


similar.  It  may  be  questioned  as  to  whether  materials  which  behave 
differently  (  Si02,  Hf02/  ZrO^)  may  be  compared  by  this  criterion. 

The  electron  diffusivity  is  given  theoretically  by  Dfi  *  ykT/e, 
where  mobility  y  *  er/m*  [Sean  1955] .  t  can  be  a  complicated  function 
of  tarperature.  This  occurs  because  it  depends  upon  the  energy  of  both 
the  electrons  and  the  phonons.  Electron  diffusivity  as  thermal 
diffusivity  is  a  difficult  quantity  to  determine  for  thin  films.  This 
quantity  is  used  primarily  to  determine  the  rate  of  spread  of  the 
electron  density.  An  order  of  magnitude  value  was  chosen  based  upon 
theoretical  considerations  for  this  order  of  magnitude  study  and  thus 
its  dependence  with  temperature  is  not  that  significant  to  the  feasi¬ 
bility  cf  the  process.  This  is  true  also  of  the  other  secondary  effect 
rate  constants. 

In  general,  tha  part  that  temperature  dependence  plays  in  the 
material  properties  is  not  that  important  to  the  proposed  process  of 
electron  density  generation.  The  temperature  effect  is  a  quantitative 
correction.  The  feasibility  of  the  process  requires  only  order  of 
magnitude  values.  The  behavior  of  the  thermal  conductivity  car.  affect 
the  quantitative  determination  of  damage  threshold.  However ,  due  to 
the  difficulty  of  a  quantitative  determination  of  the  absorption 
process  a  qualitative  (relative)  measure  of  conductivities  is  adequate. 
The  behavior  of  the  thermal  conductivity  should  net  be  important 
qualitatively.  This  is  particularly  true  if  the  behavior  is  similar 
for  various  materials.  The  outstanding  exception  of  the  similar 
behaviors  of  conductivities  is  that  of  SiO.?.  This  may,  jr.  part,  help 


explain  its  anomalously  high  damage  threshold  with  respect  to  other 
oxide  coatings. 

The  temperature  dependence  in  avalanche  ionization  may  be 
dismissed  as  well.  Avalanche  ionization  in  optical  thin  films  is  a 
secondary  process.  The  initial  absorbing  center  would  dominate  the 
initiation  of  damage  and  an  avalanche  would  simply  provide  a  mechanism 
for  the  growth  of  an  absorbing  region.  The  basic  temperature  affect  on 
avalanche  ionization  is  to  assist  in  protecting  electrons  to  higher 
energies  by  collisions  with  more  and  higher  energy  phonons. 

The  basic  conclusion  of  this  study  is  that  temperature 
dependence  of  absorption  can  play  a  part  in  the  absorption  of  optical 
radiation  in  snail  band  gap  oxide  coatings  if  avalanche  ionization  does 
-not  dominate.-  It  is,  however,  probably  a  secondary  effect  which 
results  in  the  growth  of  preexisting  absorption  center  regions. 
Temperature  dependent  absorption  cannot  play  any  significant  role  of 
absorption  in  larger  gap  coatings,  though. 


CHAPTER  8 


CONCLUSIONS 

In  Chapters  1  and  2  concepts  of  a  localized  absorbing  region 
ware  introduced.  Scaling  laws  were  derived  from  the  model  in  which 
thermal  diffusion  was  the  primary  transport  mechanise.  These  scalings 
related  the  macroscopic  thermal  material  parameters  to  the  damage 
threshold  energy  density  of  the  incident  laser  field.  The  primary 
result  from  this  work  was  that  thermal  transport  through  the  region 
surrounding  the  absorption  center  dominated  the  rate  at  which  damage 
occurred.  Damage  was  assumed  to  occur  when  sufficient  energy Jdensi ties 
were  present  to  melt  the  outside  radius  of  the  absorbing  region.  Thus, 
it  is  clear  that  the  rate  at  which  energy  diffuses  from  this  region  is 
inportant. 

In  Chapters  3  and  4,  the  generality  of  integral  transform 
techniques  was  exploited  to  investigate  other  forms  of  solutions.  In 
particular.  Chapter  3  dealt  with  the  affects  of  the  finite  boundary  due 
to  a  substrata  on  the  film  and  Chapter  4  dealt  with  repetitively  pulsed 
lasers.  It  was  showu  that  the  thermal  conductivity  of  a  substrate 
could  make  a  difference  in  damage  thresholds  of  a  film  under  certain 
circumstances .  These  circumstances  included  the  assumption  that  the 
film  and  substrate  are  in  perfect  thermal  contact.  This  may  be  a  poor 
assumption.  However,  if  the  absorption  is  at  or  sufficiently  near  the 
boundary,  electron  transport  may  occur.  Electron  transport  across  the 
boundary  is  not  the  assumed  mechanism  here.  Thermal  transport  is  the 


necbanisn  assumed.  Tbs  data  indicates  that  thermal  transport  in  the 
substrate  is  a  contributing  factor.  The  details  of  the  boundary 
condition  are  not  known.  It  was  also  shewn  that  the  laser  field 
distribution  in  the  film  can  be  a  significant  factor. 

Chapter  4  demonstrated  that  repetitively  pulsed  laser  induced 
damage  is  more  complicated  than  mere  thermal  effects.  It  also  showed 
that  thermal  build  up  in  optical  coatings  from  repeated  pulses  should 
not,  by  itself,  contribute  to  damage.  However,  a  build  up  effect 
experimentally  observed  in  optical  components  appears  to  be  due  to 
different  mechanists  in  different  materials.  The  fact  that  an 
increase,  as  opposed  to  a  decrease,  in  damage  threshold  occurs  in 
certain  films  from  the  initial  pulse  supports  this  statements'  It  is 
possible  that  photo  or  thermal  induced  changes  occur  in  the  targets. 
Thus,  subsequent  pulses  are  incident  upon  a  material  that  is  in  a 
different  state  from  its  state  preceding  the  previous  pulse.  This 
state  say  be  altered  with  every  incident  pulse  leading  to  a  build  up 
phenomena.  Such  processes  are  not  only  complicated,  but  material  ana 
structure  dependent.  A  thorough  study  of  Multiple  Pulse  Damage 
Mechanisms,  in  itself  should  prove  to  be  a  fruitful  exercise.  It  may 
very  well  contribute  to  the  understanding  of  single  pulse  damage  by 
demonstrating  acre  apparently  what  occurs  on  the  smaller  scale  of  a 
single  incident  pulse. 

Chapter  5  tied  the  findings  of  Chapters  1  through  3  to 
experimental  damage  data.  The  findings  of  these  chapters  were  shewn  to 
agree  remarkably  well  with  the  data.  Furthermore,  the  morphology  of 
damage  corresponding  to  certain  cats  \in  particular,  data  on  fluoride 


films)  was  in  support  of  the  model-  The  scaling  law  derived  in  Chapter 
2  scaled  quite  well  with  the  experimentally  measured  damage  threshold 
from  several  data  sets-  This  scaling  was  shown  to  tit  fluoride  films 
best.  A  possible  reason  for  the  poor  fit  of  oxides  is  that  the 
absorbing  region  in  sane  oxide  films  expands  to  large  dimensions 
invalidating  the  application  of  a  localized  absorbing  model.  This 
expansion  may  be  due  either  to  an  avalanche  ionization,  a  temperature 
dependent  absorption,  or  both.  This  does  not  occur  in  fluorides  most 
probably  because  of  their  larger  card  gaps.  In  addition,  experimental 
demonstration  of  increased  substrate  conductivities  providing  higher 
damage  thresholds  was  given.  This  result  was  predicted  in  Chapter  3. 

The  investigation  of  a  fundamental  mechanisn  for  providing  the 
required  absorption  was  the  topic  of  chapter  6.  Localized  absorption 
on  the  order  of  105  oa-1  does  not  exist  prior  to  the  damage  shot  and 
therefore  it  mast  evolve  during  the  shot.  A  feasibility  study  demon¬ 
strated  that  absorption  centers  provided  by  collections  of  color 
centers  were  a  dear  possibility.  Measurements  and  studies  of  absorp¬ 
tion  centers  [Schulman  and  Compton  1963]  indicate  that  colloids  may 
also  contribute  to  this  absorption  at  the  longer  wavelengths.  Multi- 
photon  calculations  of  impurity  atom  transitions  indicated  that  they 
were  an  unlikely  source  of  large  densities  of  conduction  electrons, 
though  caul ti photon  transitions  of  color  centers  nay  contribute. 

The  size  of  the  region  required  for  absorbing  sufficient  energy 
for  damage  indicates  that  it  is  likely  that  the  initial  absorbing 
region  must  expand.  Of  the  mechanises  investigated  for  this  process, 
avalanche  ionization  cf  the  off  resonant  shallow  states  (including 


inpurity  atoms)  appears  to  be  the  most  likely,  numerical  calculations 
based  upon  this  demonstrated  the  quantitative  feasibility. 

In  Chapter  7  the  effects  of  temperature  were  investigated. 

The  effects  of  both  temperature  dependent  absorption  and  the  variation 
of  material  parameters  with  temperature  were  studied.  Temperature 
dependent  absorption  was  shown  to  be  at  best  a  sectsodary  effect  for 
most  of  the  films  of  concern,  it  is  clear*  however*  that  the  tempera¬ 
ture  behavior  of  the  material  parameters  necessarily  affects  the 
quantitative  calculations.  The  qualitative  relations  should  not  be 
affected  significantly  in  most  cases  though.  Thus,  accounting  of  the 
temperature  dependent  behavior  of  the  themal  parameters  does  not  affect 
the  feasibility  of  the  processes  and  the  generalized  scalings.'  Any 
quantitative  prediction  attempts  would  require  accounting  for  tempera¬ 
ture  dependence*  particularly  when  acre  than  order  of  magnitude  values 
are  of  interest. 

What  has  been  presented  in  this  text  is  a  new  model  of  laser 
induced  damage  in  optical  coatings.  The  idea  of  an  absorbing  inclusion 
is  not  new.  The  relations  developed  from  the  concept  are*  however* 
new.  The  explicit  demonstration  of  scaling  relations  in  various  limits 
of  transport  i3  also  new.  In  addition  the  demonstration  of  the 
relationship  cf  thermal  parameters  to  experimentally  measured  damage 
has  not  previously  been  done. 

The  wavelength  dependence  observed  in  optical  coating  damage 
lead  to  the  formulation  of  a  caroot-ite  mechanism  for  laser  induced 
damage  evolution.  This  mechanism  is  also  based  upon  the  requirement  ct 
extraordinarily  high  localized  absorption.  The  mechanism  proposed  is  a 


hypothetical  construct  based  upon  meeting  requirements  of  the  findings 
from  thermal  calculations .  Cccpariscn  with  data  and  numerical  computa¬ 
tion  of  microscopic  phenomena  support  the  model's  feasibility.  It  is 
felt  that  feasibility  has  been  sufficiently  demonstrated  for  the 
process  proposed. 

The  possible  avenues  of  research  implied  in  these  findings  are 
numerous.  The  most  apparent  are:  the  difference  between  oxides  and 
fluorides,  multiple  pulse  damage,  measurement  of  the  evolution  and 
recombination  kinetics  of  generated  free  carriers,  measurement  of 
thermal  properties  of  thin  films,  the  effect  of  material  structure 
(crystalline,  amorphous,  etc.)  on  laser  damage,  measurement  of  color 
center  concentrations  in  optical  coatings,  and  the  calculation  of 
electron  phonon  and  electron  defect  collision  dynamics  in  intense 
fields,  me  answers  to  the  questions  associated  with  any  one  of  these 
areas  of  research  would  contribute  significantly  to  understanding  the 
processes  involved  in  pulsed  laser  induced  damage.  Efforts  ia  a  es^spia 
of  these  areas  of  research  is  presently  uadarwsy  at  various  institu¬ 
tions.  It  is  the  hope  of  this  author  that  this  present  work  has 
contributed  to  the  understanding  of  the  importance  of  these  various 


processes. 


APPBBIX  A 


SPHERICAL  MODEL  DERIVATION 


In  Appendix  A  the  solution  is  derived  for  the  spherical 
absorbing  region  using  the  integral  transform  method.  This  solution  is 
then  approximated  in  a  certain  region  of  behavior.  This  result  (Eq. 
(10))  has  provided  tremendous  insight  into  the  processes,  or  at  least 
the  scale  of  processes,  which  occur  in  same  cases  of  laser  induced 
damage  in  thin  films. 

As  stated  in  the  text,  the  problem  is  modeled  via  Sqs.  (Al) 
through  (AS) 
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where  toe  subscripts  i  and  h  refer  to  toe  inclusion  and  toe  host 
respectively  and  A  is  the  source  term.  The  interface  conditions 
between  toe  impurity  and  the  host  is  assumed  to  be 
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The  integral  transform  technique  involves  expansion  of  the 
solution  in  appropriate  eigenfunctions  of  the  region  which  in  this  case 
are  spherical  Bessel  functions.  Because  of  the  assured  spherical 
symmetry  we  need  only  consider  spherical  Bessel  functions  of  order 
zero.  These  may  be  written  as  sines  and  cosines  by  the  variable 
substitution  U  =  rf .  Thus  Eqs.  (Al)  and  (A2)  become 
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The  appropriate  eigenfunctions  for  this  region  are 
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where  a  and  8  are  continuous  eigenvalues. 

The  overall  tatgoral  eigenvalue  trust  be  the  sane  in  both 
regions,  i.e. 

°ia2  *  ~  *2  (A13) 

The  coefficients  Bi  and  32 ©ay  now  be  determined  from  the  interface 
conditions  for  which  we  find 
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The  ei gonfunct i oa  normalization  say  be  shown  to  be 
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which  works  out  to  be 
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with  y  2  Xa/IK  *  oa. 


The  solution  expanded  in  the  eigenfunctions  is  of  the  form 
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The  coefficients  of  the  composite  region  then  satisfy 
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where  use  has  been  made  of  the  orthogonality  of  the  composite 
eigenfunction. 

The  time  evolution  of  this  coefficient  is  found  by  applying 
this  transform  operator  to  Eqs.  (AS)  -  (A10) .  With  tne  aid  of  Green's 
theorem  the  system  and  boundary  conditions  transform  to  an  ordinary 
differential  equation  with  the  solution 
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Thus  we  find  a  solution  of  the  form 
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where  for  this  case  X  =  c  and  specifically  for  the  impurity  region 
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where  y  *  a' -  *  Sow#  if  we  assume  that  A(r,t)  *  A(t)  (i.e.-efcsorption 
within  the  region  is  independent  of  the  radius)  A  becomes 
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It  should  be  noted  that  no  assucsticns  have  been  made  as  to  the 
temporal  behavior  of  the  absorption  profile. 

The  source  term  A  is  written  in  terms  of  the  incident  intensity  I 


] 

! 


1 


in 


dr3A(r,t) 


«<¥• 


n'}I(t) 


(A30) 


where  Q  is  the  absorption  cross  section  computed  from  electrocnagentic 
theory  and  n'  is  the  imaginary  part  of  the  index  of  refraction.  The 
damage  threshold  is  then  defined  as 
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where  tp  is  the  pulse  length  required  to  reach  the  critical  temperature 
Tc  at  the  radius  a  of  the  absorbing  region.  Thus  Vp-  f(r,Q,£)  is 
inverted  to  give  E  *  g (T^,tp,r,Q)  for  the  damage  threshold  in  Joules 
per  centimeter  squared. 

Now,  for  a  single  square  .modeled  pulse  the  integral  becomes 
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The  first  teen  is  shown  to  be  analytically  soluble  by  the  Laplace 
transform  approach  to  this  problem  (Goidenberg,  Tranter,  1952)  or  by 
multiple  integration  by  parts  when  c  *  0  and  b  *  1. 

Thus,  assuming  Q  *  ^a2 land  E  =  Itp, 
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Since  for  the  region  'where  t ^/y  >  1  the  exponential  dominates.  The 

■y 

F{y)/y"  may  be  expanded  in  a  Taylor  series  about  y  *  0  and  integrated 
term  by  term.  The  first  term  of  the  Taylor  series  is,  for  example, 
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As  stated  in  the  text,  the  damage  threshold  has  a  minimum 
versus  radius.  That  is,  there  is  a  size  of  absorbing  region  most 
likely  to  dasage.  This  is  found  by  taking  (dE/da)  =0.  In  the  first 
order  approximation  above,  this  is  found  to  be 
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Expanding  to  the  next  order  and  integrating,  the  solution  is  found  to  be 
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where  the  second  term  in  braces  is  usually  quite  srall.  Thus  we  nay 
write 
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Equation  (A37)  may  be  obtained  directly  by  skipping  the  next  order 
integration  which  lead  to  Eq.  (A36) .  Performing  the  next,  order 
integration  allows  one  to  see  trie  validity  of  neglecting  it. 


APPEIOIX  B 


CYLINDRICAL  MODEL  DERIVATION 

In  Appendix  B  the  solution  is  derived  for  the  cylindrical 
absorbing  region  using  the  integral  transform  method.  This  solution  is 
then  specified  to  a  X/2  standing  wave  between  an  insulating  surface  and 
an  insulating  substrate. 

Mathematically  the  model  is  described  by  the  following  thermal 
diffusion  equations: 
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where  A  is  the  source  term  in  watts  per  cubic  centimeter. 

Because  there  are  no  expeincental  data  with  which  to  assess  the 
thermal  impedance  between  the  host  and  the  inclusion  (i.e.  the  heat 
transfer  coefficient),  ideal  boundary  conditions  are  assumed,  i.e. 
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Ki  dz  KS  dz 


at0^r<a,z=i 


T.  =*  T 
u  s 


dTh  dTs 

:<h  ar  *  Rs  ar  at  a  4  £’a  * 1 


at  2  *  0 


Poe  the  special  cases  of  the  infinite  and  the  nonconducting 
substrates  the  prcblen  reduces  to  Eq.  (31)  and  (B2)  together  with  Eq. 
(B5)  where  Eq.  (B6)  and  (37)  gotoT=»0atz=»l  for  an  infinitely 
conducting  substrate  and  cT/dz  *  0  at  z  *  1  for  a  nonconducting 
substrate.  Exact  solutions  are  found  for  these  cases. 

It  should  be  noted  that  cany  cccaacn  substrates  employed  in  high 
energy  laser  reflectors  are  very  poor  conductors  (e.g.  Si02) .  Further¬ 
more,  as  was  shown  for  ti>e  carcsn  half-wave  film  the  temperature  at  the 
boundary  with  the  substrata  does  not  change  significantly  for 
sufficiently  short  pulse  lengths.  The  case  of  the  nonconducting 
substrate  therefore  becomes  an  important  one  to  consider. 


The  solution  to  this  particular  case  may  be  developed  by 
considering  corresponding  homogeneous  equations  for  the  combined 


system.  The  obvious  eigenfunctions  for  the  regions  are 

v-  »  J0(or)  cos(S|S)  (99) 


and 


^  =>  lAJ0(8r)  +  3Y0(sr)}  cos(^)  IBiO) 

where  J  and  Y_  are  the  Bessel  functions  of  the  first  and  second  kind 

0  v 

respectively  and  of  order  zero,  and  a  and  8  are  eigenvalues.  Thus,  we 
find  solutions  of  the  fora 
_  1 

TR(r,t)  *  M  do  a  cnta,t,*Wa'r)  R  *  i,h  {Bll) 

n  J 

where  the  overall  temporal  eigenvalue  oust  be  the  same  in  both  regions, 
i.e. 


Di{a2  +  fp)2}  »  D.q{82  -  (y~) 2}  (312) 

The  coefficients  A  and  8  may  easily  be  solved  for  by  satisfying 
bound?3ry  condition  (B6)  for  the  functions  and  their  derivatives.  We 
find  that 


A  * 


ti^aJ^aaJYglSa)  -  ^SJq  (oa^  (3a) } 


(3125 


and  the  eigenfunction  normalization  may  be  shewn  to  be 


U  (a)  =  ^  +  Bn2(a>*  ?T  <315) 

n  i  n  n 

The  coefficients  (^(a tc)  may  be  found  via  orthogonality  as 

/  rl  k— 

CL(a.t)  *  y4rr  >  f  dr  r  f  dz  3  *  (or)T  (r,t).  (316) 
a  *Va'  R*i,h  it  'G  °R  80  R 

Thus,  to  solve  for  Cn(a,t)  explicitly  the  above  transform  is  applied  to 
the  governing  differential  Eqs.  (Bi)  and  (B2) .  With  the  aid  of  Green’s 
theorem  the  systar  and  boundary  conditions  transform  to  an  ordinary 
differential  equation  with  a  solution  of 

a  2  nrr  2  (a,t) J, (aa) 

Co",'tl  *  “S’HMc2  Hf>  )tl  — -  (Bi7) 

where 


P dr*  f  dr  cos(%  exp  ID  -  (a^  +  (BI)2}t’]A(z,t»)  (B18) 
'0  J0  A 


The  resulting  solution  for  the  temperature  at  some  point  in  the 
impurity  and  at  a  certain  time  then  becomes 


0 

T-  (r,t)  »  I  !  dx  a  C  (a,c)J  (ar)COS(£p)  (BIS) 

n  ' 

The  source  term  A(z,t’ )  which  describes  the  spatial  and  temporal 
absorption  of  radiation  within  tha  inclusion  is  quite  general,  ijo 
functional  dependence  on  z  and  t  has  been  assured.  This  solution, 
therefore,  can  take  into  account  the  spatial  dependence  of  standing 
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wave s  in  the  film  and  sue*  temporal  dependence  as  a  repetively  pulsed 
laser.  Any  spatial  dependence  can  be  repetitively  pulsed  by  forming  a 
sue  of  integrals  over  time,  integrating  fran  (n  -  l)t  +  (n  -  l)t,  to 
nt^  +  (n  -  l)td  v^sere  t^  is  each  pulse  length,  t^  is  the  delay  time 
between  pulses  and  n  is  the  pulse  nusber. 

For  the  realistic  case  of  a  single-pulse  sine- squared  standing 
wave  (half-wave  film)  the  source  term  is 


A(z,t)  =*  Aq  sin2  (-p)  *  ~  {1  -  cos {—) } 
In  this  case  the  coefficients  became 


(B2Q) 


Cn(a,CJ  ~  2H  (a) 


Ag3  1  -  e3q>{-Di(a^+(ns/l)2}t} 


Dila2  +  (nsA)2} 


2_  ^(aa) 

'5n0  "  5n2} 


(321) 


so  that 


T-  *  Y  ^  .  *  T--  +  T_-  (B22) 

i  L  m  Oi  2i  v  ' 

n 

The  source  term  A(z,t)  can  also  be  written  in  terms  of  the 
incident  laser  intensity  I (W  cm-2)  and  the  absorption  cross  section 
Q(oa2)  by  realizing  that 


dr"  A(r,t)  *  \il  (t)  0  <  t  <  t 

!  _  D 

*  R 


A  »  0 


t  <  t 


(323) 


wnere  t  is  the  duration  of  tne  Laser  ouise-  In  this  case 


~  -a2l  -  IQ 


or 
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(B24a) 


(B24b) 


&S5?£?©IX  C 


REPETITIVE  PULSE  DERIVATION 

In  Appendix  C  the  process  for  extending  all  solutions  in  this 
text  to  the  multiple  pulse  case  is  demonstrated.  By  integral  trans¬ 
form,  one  is  left  with  an  ordinary  differential  equation  in  time  wnich 
has-  the  general  solution  of 


A  »  erp(-Ct)  f  dt* 
JQ 


(Cl) 


(see,  for  example,  Eq.  (39) ) .  For  the  repetitively  pulsed  case  where 
only  thermal  effects  result. 


A  »  exp(-^t)F(a)A  {f^dt*  +f 

Jo  Jt 


t,+  2t 
d  p 


dt*  e4t’ 


“d+tp 


*  ... 


(S-l) (td+tp)  ♦  t£ 
(N-i) (td+tp) 


dt’  e$t’} 
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t  -  +  t  *  - 
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*■  J? *1 

tA 

i*I 


multiply  by  e 


fi-~k  T  (eCT)n 

*  n»Q 


(e  p  -  1)  (1  -  e^T) 
C  (1  -  e"T) 


_  ^[(N  ^  1)T  +  V 


to  find 


—  C-  -  e  {1  -  e~^T) 

A  *  (A  F(a) )  (-^ — | - L  — — —r2-) 
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